UNCLASSIFIED 

405  153 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


40  5  X  5  3 


AFSWC-TDR-63-2 


0-  /_  tr 


swc 

TDR 

63-2 


STEP  LOAD  MOVING  WITH  LOW  SUBSEISMIC  VELOCITY 
ON  THE  SURFACE  OF  A  HALF- SPACE  OF  GRANULAR  MATERIAL 


April  1963 


TECHNICAL  DOCUMENTARY  REPORT  NUMBER  AFSWC-TDR-63-2 


\ 


Research  Directorate 

AIR  FORCE  SPECIAL  WEAPONS  CENTER 
Air  Force  Systems  Command 
Klrtland  Air  Force  Base 
New  Mexico 


Project  No.  1080,  Task  No.  108001  ^ 

U  T1SIA  & 


This  research  has  been  funded  by  the 
Defense  Atomic  Support  Agency  under  WEB  No.  13. 145 


1 


(Prepared  under  Contract  AF  29(601)-2855 
by  Hans  H.  Bleich  and  Ewald  Heer,  Paul 
Weidlinger,  Consulting  Engineer,  New 
York,  New  York) 


i 


HEADQUARTERS 

AIR  FORCE  SPECIAL  WEAPONS  CENTER 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the  Government 
may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation, 
or  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  is  made  available  for  study  upon  the  understanding  that  the 
Government's  proprietary  interests  in  and  relating  thereto  shall  not  be  im¬ 
paired.  In  case  of  apparent  conflict  between  the  Government's  proprietary 
interests  and  those  of  others,  notify  the  Staff  Judge  Advocate,  Air  Force 
Systems  Command,  Andrews  AF  Base,  Washington  25,  DC. 


This  report  is  published  for  the  exchange  and jstimilation  of  ideas;  it  does 
not  necessarily  express  the  intent  or  policy  of  any  higher  headquarters. 

Qualified  requesters  may  obtain  copies  of  this  report  from  ASTIA. 

Orders  will  be  expedited  if  placed  through  the  librarian  or  other  staff  member 
designated  to  request  and  receive  documents  from  ASTIA. 


TDR-63-2 


ABSTRACT 

The  plane  strain  problem  of  a  step  load,  p,  moving  on  the  surface  of  a 
half- space  of  granular  material  governed  by  internal  friction  and  cohesion  is 
being  considered  because  of  its  importance  to  Air  Force  research  in  pro¬ 
tective  construction.  As  a  preliminary  step  towards  a  more  complete  solution, 
in  this  paper  it  is  assumed  that  in  regions  of  slip  the  elastic  deformations  may 
be  neglected  in  favor  of  those  resulting  from  slip.  The  latter  assumption 
limits  the  application  to  Mach  numbers  of  less  than  about  0.  20. 

Two  possible  types  of  behavior  of  the  material  are  considered.  During 
slip,  one  material  exhibits  dilatancy,  while  the  other  does  not  change  in 
volume.  Because  of  neglect  of  the  elastic  deformations  in  the  slip  region, 
only  degenerate  results  are  obtained  for  the  case  of  dilatancy;  while  stresses 
can  be  determined,  deformations,  velocities  and  accelerations  vanish.  For 
the  other  material,  all  desired  quantities  are  obtained. 

A  significant  finding  is  that  for  values  of  p  above  a  certain  limit,  granular 
particles  will  be  expelled  at  the  surface  ahead  of  the  pressure  front.  The 
applied  pulse  will  be  preceded  by  a  precursor  of  expelled  grains. 
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Throughout  this  report,  equations  are  numbered  by  hyphenated  numbers,  such 
as  Eq.  (2-7),  (3*18),  or  (B-lb).  The  first  number,  2,  3,  B,  respectively. 
Indicates  that  the  equation  occurs  in  Section  2,  3,  or  in  Appendix  B. 

When  an  equation  is  referred  to,  the  full  number,  say  (2-7),  Is  quoted  If 
the  reference  occurs  in  a  different  section.  However,  if  an  equation  is 
mentioned  in  the  section  in  which  it  originally  appears,  only  the  seoond 
half  of  the  number  is  quoted 1  In  Section  2,  Eq.  (2-7)  would,  therefore, 
be  quoted  simply  as  Eq.  (7),  etc. 
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1.  INTRODUCTION 


Th«  present  paper  la  tha  first  step  toward  the  determination  of  tha 
•ffaeta  of  a  praaaurt  distribution  p  (x-Vt)  moving  with  a  velocity  V  on  tha 
surface  of  a  half-space  consisting  of  a  granular  material,  (Pig.  1),  Because 
of  the  complexity  of  tha  matter  only  the  steady-state  situation  Is  being 
studied,  1.  e,|  the  effects  of  Initial  (starting)  conditions  are  excluded. 

The  equivalent  steady-state  problem  for  an  elastic  half-space  has  been 
treated  by  Cole  sad  Huth,  [Ref.  1),  for  the  case  of  a  moving  line  load 
p  (x-Vt)  •  4  (x-Vt),  where  4  Is  Dirac's  delta  function.  Using  the  principle 
of  superposition,  the  effect  of  an  arbitrary  load  p  (x-Vt)  can  be  found  by 
integrations.  Considering  the  behavior  of  a  granular  material  to  be  governed 
by  internal  frlatlon  of  the  Coulomb  type  and  possibly  cohesion,  the  elastle 
solution  can  be  applied  up  to  a  certain  intensity  of  the  loada  pt  beyond  this 
critical  intensity,  which  depends  on  the  velocity  V  and  the  distribution  of  p. 
Internal  slip  must  occur  In  the  material  requiring  a  different  analysis  which 
Is  to  be  developed  la  this  paper. 

Por  a  complete  understanding  of  all  possible  situations,  the  elastle 
properties  of  the  material  and  the  slipping  effeots  must  be  taken  Into  account 
simultaneously.  However,  the  problem  becomes  very  Involved,  and  It  appears 
advantageous  to  treat  at  first  the  simple  esse  of  a  stop  loaA,(Plg.  2) ,for  an 
elastically  extremely  rigid  material,  such  that  the  elastle  deformations  can  be 
assumed  to  vanish.  In  spite  of  the  serious  Implications  of  this  assumption, 
which  limits  the  application  of  the  theory  to  low  Mach  numbers,  the  results 
obtained  permit  a  physical  understanding  of  a  peculiarity  of  the  Coulomb  type 
material  in  dynamic  situations i  When  the  pressure  wave  exceeds  a  certain 
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intensity,  a  pm cur* or  occur* , con* 1* ting  of  particle*  expelled  from  tha 
surface  ahead  of  the  applied  pressure  wave. 

The  equivalent  of  the  dynamic  problem  treated  here  Is  the  static  one 
of  determining  the  limiting  load  p  which  may  act  on  one  half  of  the  surface 
of  a  half-space,  (Fig.  3a).  This  problem  la  a  special  ease  of  one  treated 
by  Prandtl,  (Ref.  2],  who  considers  a  load  p  of  finite  length  1,  (Fig.  3b). 

The  present  treatment  for  the  limiting  case  V  ♦  o  overlaps  Prandtl's  problem  for 
1  - 


In  order  to  formulate  equations  of  motion,  the  "yield"  condition  between 
the  stresses  defining  tha  occurrence  of  slip  must  be  supplemented  by  a  flow 
rule  defining  the  strain  rates.  Assuming  the  material  to  be  Isotropic,  the 
principal  axes  of  the  stress  tensor  and  the  strain-rate  tensor  must  necessarily 
coincide  everywhere.  However,  this  statement  dose  not  sufflee  end  further 
relations  are  required.  In  accordance  with  the  theory  of  plastlelty,  Drueker 
and  Prager,  (Ref.  3],  have  formulated  such  relations  by  postulating  the 
existence  of  a  plastic  potential.  The  condition 
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In  terms  of  the  invariants  j»  Jj  *  J  ((  •  eg)8  ♦  (eg  -e^8 

is  the  generalisation  of  the  Coulomb  yield  condition  for  a  three-dimensional  state 

of  stress,  where  o^,  o^,  a ^  are  the  principal  stresses  and  where  s  end  e  are 

appropriate  constants  connected  with  the  values  of  the  slip  angle  and  cohesion, 

[Ref,  3].  The  concept  of  the  plastic  potential  then  defines  the  plastic  strain 

rates  I  as  derivatives  of  the  yield  function  f  with  respect  to  the  stresses 
•J 

V 


• *  if 


(1-2) 
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where  1  Is  an  arbitrary  poaitiva  quantity  vhieh  nay  ba  a  function  of  tba 
location,  for  tba  eaaa  of  plana  atrainlof  intaraat  here,  tha  resulting 
stress-strain  rata  ralationa  ara  given  in  (Ref.  •»],  It  la  noted  that  tba 
uaa  of  tha  yiald  condition  f  aa  plaatle  potential  rcqulrca  that  clip  ba 
accompanied  by  a  rolune  exp ana loo. 


A a  aa  alternative  to  tha  ideal  iaotrople  material  defined  by  [}],  one 
can  define  tha  behavior  of  a  material  in  plana  atrain  completely  by  supplementing 
the  Coulomb  yiald  condition ,and  tha  requirement  that  tha  ana  of  tha  etreee 
tensor  and  atrain  rata  tensors  agree ,by  tha  statement  that  tha  plaatle 
deformations  ara  lncoapresalbla.  Thie  aaaumptlon  has  baan  utlllsad  by 
Iahllnaki ,[Raf.  5], to  formulate  aquations  of  motion  for  plana  problems  and 
for  problems  of  spherical  symmetry.  Tha  assumption  of  incompressibility  which 
appears  not  unreasonable  -  or  physically  impossible  -  is  incompatible  with 
the  formulation  by  Drunker  and  Pragar.  To  clarify  tha  fundamental  difference 
between  the  two  materials,  one  can  generalise  the  "assumed"  Incompressibility 
in  plane  strain  to  a  three  dimensional  statement,  for  an  Isotropic  material 
any  condition  nay  be  only  la  tame  of  Invariants |  tha  appropriate  flow  rule 
defines  the  strain  rates  again  as  derivatives  of  a  plastic  potential,  but 
tha  function  to  ba  used  is  f^  ■  F  (Jg),  where  f  may  ba  an  arbitrary  function. 

Aa  tha  result  doea  not  depend  on  tha  form  of  the  function  F,  its  choice  only 
affaots  tha  value  of  the  arbitrary  factor  1,  one  may,  without  loss  of  generality, 
uaa  f j  •  Jg.  Tha  combination  of  lq.  (1)  and 


•  X 


M. 

T 

44 


(1-3) 


defines  tba  three-dimensional  behavior  of  the  material  postulated  la  [?]• 
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The  difference  in  the  formulations  of  [3]  and  [$]  lies  In  the  fact  that 
(3]  assumed  tba  function  f  in  th#  yi#ld  condition  to  ba  th#  plastic  potential. 

It  1#  conventional  to  us#  th#  yield  function  as  plastic  potential  but  th#  us# 
of  another  function  of  th#  invariants  will  not  b#  ruled  out  here,  in  spit# 
of  certain  thermodynaaic  difficulties  resulting  froa  th#  us#  of  different 
functions*)  in  th#  yield  condition  and  in  th#  flew  rul#. 

Without  attempting  to  decide  on  th#  physical  applicability  of  the  two 
concepts,  th#  following  paper  will  pursue  both  possibilities.  It  will  be  seen 
that  th#  solutions  of  dynaaieal  problems  for  th#  dilating  material  (Drucker-Prager) 
become  degenerate  if  the  elastic  deformations  ar#  Ignored,  such  that  the  alters*- 
tlv#  case  of  the  "incompressible"  material  gives  more  interesting  result#  and 
its  analysis  will  be  emphasised. 

A  comment  on  steady-state  solutions  of  the  type  considered  here  is  In 
order.  Due  to  the  fact  that  the  manner  in  which  the  loads  at  large  distances  from 
the  front  have  been  applied  at  early  times  is  not  Included  in  the  statement  of  the 
problem,  there  can  be  no  uniqueness  or  existence  theorems,  and  there  m y  be  BK>re  them 
one  solution  for  a  given  case.  In  the  elastic  problem  [1],  there  is  no  steady-state 
solution  if  the  velocity  V  equals  that  of  Rayleigh  waves}  for  other  velocities 

•j  An  elementary  consequence  of  the  two  different  assumptions  may  be  of  interest. 

In  the  non-dilating  material,  Eq.  (3),  in  the  csss  of  plane  strain  the  transverse 
stress  0^  is  equal  to  the  mean  of  the  other  two  principal  stresses 

3  2 

For  the  dilating  material,  on  the  other  hand,  the  absolute  value  of  the  trans¬ 
verse  stress  is  necessarily  larger  than  the  mean, 


the  excess  depending  on  the  state  of  stress  and  on  the  material  constants, 
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the  solutions  still  contain  arbitrary  constants.  The  horlaontal  strata 
contains!  e.g.,  an  opan  constant,  equivalent  to  tbs  superposition  of  a  uniform 
stats  of  horlaontal  stress)  this  lndeteralnaoy  can  not  he  removed.  In  some 
eases,  solutions  can  be  excluded  by  additional  general  considerations  of  the 
starting  situation  one  has  in  mind.  This  is  best  seen  on  the  elementary  example 
of  a  half-space  of  an  inviscid  compressible  fluid,  loaded  by  a  pressure  pulse  p 
progressing  vltb  supersonic  velocity  V,  There  is  an  obvious  solution,  (Fig,  ha), 
where  the  load  produces  a  plane  wave  of  intensity  p,  progressing  vltb  a  front 
inclined  at  an  appropriate  angle  a.  However,  this  le  not  the  only  steady-state 
solution)  an  alternative  is  a  plane  wave  whose  front  is  inclined  at  an  angle 
of  180  -  e.  Combinations  of  the  two  are  also  eorreot  steady-state  solutions. 

If  one  is  interested  in  states  generated  by  the  application  of  pressures  on  the 
surface  only,  one  can  reason  that  solutions  Including  the  wave  front  shown  In 
(Fig,  bb)  can  not  occur,  such  that  a  unique  solution  Is  obtained. 

Similar  simple  considerations  do  not  exist  in  subsonic  situations,  leaving 
one  even  In  the  relatively  simple  elastic  case  with  non-unique  solutions.  It  is, 
therefore,  to  be  expected  that  non-uniqueness  may  also  occur  in  the  case  of 
the  material  considered  la  this  paper. 


5 


2.  BEHAVIOR  OF  MATERIAL 


It  will  be  found  convenient  to  formulate  tb«  equations  describing  tba 
behavior  of  tba  material  In  an  unuaual  manner,  using  tba  major  principal 
stress  and  its  direction  as  independent  variables. 

In  accordance  vith  tbe  usual  conventions ,  positive  stresses  and  positive 
strain  rates  indicate  tension  and  elongation,  respectively.  The  inequalities 
required  to  describe  tbe  material  behavior  will  be  on  tbe  quantities  -s^, 

-Oj,  eto.  As  tbe  material  to  be  considered  requires  that  at  least  tbe  major 
principal  stress  *|  is  compression,  it  follows  that  tbe  quantities  at  each 
side  of  tbe  inequalities  are  positive  numbers. 

Two  types  of  material  will  be  considered.  Both  are  assumed  to  permit 
no  elastic  deformations  at  all|  any  deformations  vhlob  occur  are  caused  by 
slip.  Whenever  such  slip  occurs,  tbe  Coulomb  condition  between  tbe  streaees 
must  be  satisfied. 

CASE  A.  Incompressible  Material f 

This  material  agrees  with  tbe  one  assumed  in  ($)  and  is  described  by 
tbe  following  assumptions  I 

1.  The  material  is  incompressible,  i.  e.,  in  addition  to  assuming  that 
tbe  elastio  deformations  vanish,  there  is  no  volume  change  during 
slip. 

2.  Tbe  state  of  stress  wbicb  si/  exist  in  any  element  is  restricted, 
such  that  at  least  ons  of  tbe  two  principal  stresses,  (Fig.  ?), 
must  be  compressive |  1st 

-oj  >  -o^  (2-1) 
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Mid  the  restriction  requires 


•a,  >  Oj  (2-2) 

9}  viU  be  referred  to  aa  tba  ‘‘major"  principal  atraaa. 

J.  Tba  permissible  atataa  of  stress  ara  further  restricted  by  the 
condition i 

•*l  i  ko  (2-3) 

vhere  k  la  i  poaltire  number,  k  <  1,  which  defines  the  ancle  of 
interior  friction  and  the  variable  0  la  deflced  by 

0  ■  -(0j  ♦  a)  (2U») 

In  this  relation  a  >  0  la  a  materiel  constant  defining  the  amount 
of  cohesion.  The  reguiresMnt  0,  «  0  leads  to  the  eondition  o  >  -at 
if,  on  physical  grounds,  one  wishes  to  Insure  that  ot  is  also 
compressive,  0  >  0  could  be  specified, 
k.  If 


•a I  *  ko  (2-5) 

the  eleaent  will  be  rigid,  and  all  strains  vanish.  If,  however, 

•a2  *  k»  (2-6) 

the  element,  while  incompressible,  may  deform  but  only  in  such 
fashion  that  the  principal  axes  of  the  strain  rate  tensor  t|,  s2, 
are  parallel  to  o,,  o2,  respectively.  In  addition,  the  strain 
rate  c{  must  be  negative  or  vanlshi 

«,  *  0  (2-7) 

Due  to  the  assumed  incompressibility 

*,  ♦  c2  ■  0  (2-8) 

such  that  Ig.  (7)*  assures  dissipation  of  energy  during  deformation. 


IT  fees 'note  on  page  ii. 
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The  ititnuti  1.  to  >t.  express  the  situation  la  a  material  with 

a  coefficient  f  of  lntarnal  Coulomb  friction  and  cohesion.  Tbs  relation 

batvaan  k  and  f  can  ha  found  in  tha  following  manneri  Consider  a 

rectangular  element,  (Fig.  5) ,  under  tha  action  of  tha  principal  stresses 

Oj  and  o2.  If  a  plane  Baking  an  arbitrary  angle  a  with  the  direction  of 

o.  is  drawn,  the  normal  and  tangential  stresses,  a  and  t  •  respectively, 
1  so 

are  t 

o  ■  o  cos*a  ♦  0  sin*o 
a  1  2 

t  ■  (0  •a  )ain  a  cos  a  (2-9) 

a  i  2 


Mo  slip  can  occur  if  for  all  values  of  a  the  direct  pressure 

Bultlplled  by  the  friction  coefficient,  -fo  ,1s  larger  than  the  excess  of 

a 

the  shear  stress  |i  |  over  a  value  |i  |  which  defines  cobsslon,  or 


!'•!  -  lT-l 


<  f 


The  condition  of  slip  will  be  reached  if  the  largest  value  of 
the  above  ratio  as  function  of  a  just  reaches  f,  l.e.,  when 


l'°l  , 


(*-M) 


Substitution  of  0^  and  from  Eqs.  (it)  and  (6)  into  Eqa.  (9) 
gives  after  rearrangenent 


o  »  -e  ~  -o  tLzJl  cos  2a  *  s  coaza 


JTa|  ■  c  -  2  ~  *1®  2a  ♦  jf  2a 


(2-11) 
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To  find  the  anglo  a  for  which  the  naxinun  of  the  fraction  In 
Kq.  (10)  occur*,  Kqs.  (11)  arc  substituted,  and  differentiation  with 
respect  to  a  is  performed.  After  cancellation  of  canon  factors,  one 
finds  a  condition  for  at 

[o  ”2  ^  ♦  «  ■  fe""  cos  2e  ♦  ■  cos^aj  cos  2a  ♦ 

♦  |o  sin  2a  ♦  |  sin  2a  -  | ’s |J  sin  2a  •  0  (2-1 2) 

Noting  from  Kq.  ( 11)  that  the  expressions  in  parentheses  are, 
respectively,  the  nunsrator  and  denoolnator  of  Kq.  (10),  one  obtains  f  In 
terms  of  the  ns  yet  unknown  angle  si 

fa. 

•in  2a  (2-13) 

Rearranging  Sq.  (12)  to  separate  terns  which  contain  0'  and  those 
which  do  not,  gives  after  simplification 

o g  -  00s  2s  ♦  l—g ~j ♦  cos  e  [.  cos  o  •  8|fg|  sin  ej  •  0 

This  relation  nust  hold  regardless  of  the  value  of  s,  requiring 
that  the  coefficient  of  0  and  the  second  tern  each  vanish  individually. 
The  vanishing  of  tbs  coefficient  of  0  gives  a  condition  on  si 

CM  .  $-=-£  (2.15) 

from  Kq.  (13) 

f  .  .£SSJ£_  mLsjL 

/l  -  cos*2s  2rit  (2-16) 

The  vanishing  of  the  second  tern  in  Kq.  (Ill)  relates  finally 


9 


| Ts t  and  at 


IM 


•  a 


2  1  ♦  eoa  2a 

vhere  the  identity  cot  a  ■  ^  Z  Coa  2a  ***  UM** 


(2-17) 


CASE  B.  Material  vith  Dilat&ncy 

It  ia  again  aaauaed  that  tha  alaatie  daforaationa  of  the  Material 
vanlah,  but  in  accordance  vith  tha  aaaumptlona  in  (3]  there  vill  be  an 
increaae  in  volune  aaaociatad  vith  alip. 


Retaining  Eqa.  (3)  and  (b),  tha  naterial  vill  reaaln  rigid  if  tha 

inequality  (5)  applied,  and  nay  alip  if  Eq.  (6)  ia  aatiafied.  In  thla 

case,  the  direction  of  tha  principal  atrain  rataa  c(  and  eolnclda  vith 

thoae  of  the  principal  atreaaaa  oJ#  raapeetivaly.  Tha  valuaa  of  i,, 

c,  required  by  the  aaaunptlona  in  (3)  for  tha  eaaa  of  plana  atrain  can  be 

found  conveniently  froa  (Ref.  k,  Sq.  6)  by  letting  tha  a,  y-axea  coincide 

vith  the  principal  anea,  and  noting  that  aln  i  in  tha  reference  haa  the 
1  -  k 

value  j  in  tha  ayabola  employed  herai 

J, -KHr*1]  «*»> 


vbara,  alnilarly  to  Sq.  (1-2),  1  danotaa  an  arbitrary,  location  dependent, 
poaltlve  quantity. 
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3.  SOLUTIONS  IN  REQIOMS  WITH.  AMD  WITHOUT  SLIP.  RESPECTIVELY . 


As  a  preliminary  to  tbs  construction  of  solutions  satisfying  tbs  boundary 
conditions  on  tbs  surfaes,  it  Is  convenient  to  consider,  respectively,  solutions 
in  regions  vbsrs  sitbsr  Kq.  (2-5)  or  Xq.  (2-6)  Is  satisfied,  l.e.,  rsglons 
without  slip,  and  with  slip,  respectively. 

In  areas  vbsrs  Iq.  (2-6)  is  satisfied,  tbs  squatlons  of  notion  in  Cartsslan 
coordinates  x,  y  with  rsspsot  to  tbs  fixed  origin  0,  (fig.  6),arsi 


Jo  It 

Ju  Ju  Jul 

-  ♦  -  a  p 

— -  ♦  u  —  ♦  v  — 

ax  ay 

at  ax  lyl 

Jt  lo 

Jv  It  9vl 

- ♦  a  g 

—  ♦  u  —  ♦  v  — 

ax  by 

L»*  i*  iyj 

vbsrs  u,  v  su«,  rsspsetlvely ,  tbs  x  and  y  oonponsnts  of  tbs  velocity,  and  e  ,  • 

—  w 

and  T  arc  tbs  strsssss. 


Consider  an  element,  (fig.  5),  in  vhleh  tbs  major  principal  stress  is 
inclined  at  an  angle  0  to  tbs  horisontal*  In  tbs  condition  of  slip,  Bqs.  (2-b) 
and  (2-6)  give 

e)  a  -o  -  s 

o,  ■  -  be  ( 3-2) 

a 


vbsrs  b  and  s  are  positive  quantities,  expressing  tbs  stress  sonponents  in  tens 
of  tbs  prinolpal  stresses  and  of  the  angle  I,  one  obtains  after  substitution  of 
(2) 


0 

X 


s 

y 


t 


+  cos  2eJ  -s  cos*# 

-  idl  cos  2lJ  -e  sin*# 
sin  21  -  £  sin  29 


(S*3) 


li 


From  this  point  on,  tbs  tvo  types  of  mats rial  dasoribsd  in  Ssetion  2 
require  ssparats  trsatasnt. 


A,  Region  with  Slip  in  the  Incompressible  Material 

As  stated  in  Ssotion  2,  the  principal  strain  rates  and  t2  are  parallel 
to  Oj  and  o2f  further,  froa  Iq.  (2-7),  <}  *  0.  Introducing  a  new  variable, 
the  positive  quantity  a ,  defined  by 

e  i  (3-M 

the  condition  (2-8),  expressing  incompressibility,  requires 


«  ■  e 


(3-5) 


The  components  of  the  strain  rate  tensor,  i  ,  c  ,  c  ,  can  now  be 

**  v  /  "v 

expressed,  alternatively,  la  terns  of  §  and  •,  or  in  terms  of  u  and  v*  This 
leads  to  the  following  relations i 


*w  (3-6a) 

Jyy  ■  asa  21  ■  ^  (3-6b) 


At  this  point,  use  is  made  of  dimensional  considerations  and  of  the  fast 
that  only  steady-state  solutions  are  to  be  determined.  Stresses,  velocities, 
strain  rates,  etc.,  can  be  expressed  in  the  general  form 


yna 


'  E5*  f  M 


where  f  means  a  function  of  the  variables  stated. 


(3-7) 


If  the  ratios  &  and  -JU  are  considered  as  parameters  defining  the 

s  py* 
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particular  physical  problem,  the  exponents  are  uniquely  defined  by  the 
dimension  of  the  expressed  quantity.  Therefore i 


o.  * 


•y  »  T  •  »  ' 


p  f(l) 


?u,  Vv  ...  -  J  f(f) 


Ve 


■  ;  r  ,(t> 

•  f(f) 


(3*8) 


where 


( 


x-Vt 

y 


(3-8a) 


and  the  functions  f  differ,  of  course,  for  different  quantities.  The  derivatives 
of  f  ■  f ( y-)  I  f(|)  with  respect  to  x,  y  and  t  are 


JMT  1  df 
»x  y  df 


If  .1  df 

F  y  « 


(3-9) 


Applying  these  relations  to  the  derivatives 


in  Iqs.  (1)  and  (6)  gives 


§■-«$— HI 


(3-10) 
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and 


du  • 

■  -  ye  eea  21 

C  •  -  ye  cos  2d  (3-11) 

dv  _  du  _  •  . 

dT  -  <  ft  •  -  •*»  2» 

Ttaa  tarm  a  balng  a  function  of  C  which  night  vanish  Identically ,  tha 
tvo  aituatlona  a  (I  0  and  a  ■  0  require  separata  canal  deration* 

Case  1.  a  4  0 

Bllmlnation  of  ^  and  ^  in  Eqs.  (11)  furniahaa  tha  relation! 

( C  -  i)  eoa  21  ♦  2  ala  21  ■  o  (3-12) 

The  variable  (t  defined  by  Kq,  (8a) .can  ha  expressed  by  the  angle  0, 

(Mg.  6).  defining  the  position  of  the  elensnt  with  respect  te  the  noving  front! 

<  •  aot  0  (3-13) 

Substitution  into  Kq.  (12)  gives  tvo  possible  values  for  the  direction 
of  the  najor  principal  axis 


•  ■  0  ♦  f 


(3-l4) 


To  obtain  the  stresses,  one  finds  froi  the  first  two  Bjo.  (11) 

which  pernits  elimination  of  all  velocity  terns  from  Kqs.  (10), 


(3-13). 


Changing  the  independent  variable  la  the  resulting  differential  equations 
from  (  to  0,  one  obtains 


do*  d»  a 

fe*  -  9  ♦  ®«t*0 


5&1  & 

d0  I  d{ 
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(3-16) 


Excluding  the  points  0  ■  0,  »  vhtrs  ^  vanishes,  substitution  of  Eqs.  (3) 
furnishes  a  first  order  differential  equation  for  the  quantity  ot 

(0  0  0,  •)  (3-17) 

Its  solution  is 

(J.l8) 

The  upper  or  lover  signs  are  to  be  selected  In  the  sene  meaner  as  in 
Eq.  (lb), while  C  is  an  arbitrary  constant. 


The  accelerations  a^  and  ay  ,  and  velocities  u  and  v  can  be  obtained  from 
Eqs.(l)  and  (10).  The  right  hand  sides  of  Eqs.  (l)  are  pa^  and  pa^,  respectively, 
and  one  recognises  therefore  that  the  right  hand  sides  of  Eqs.  (10)  are  pya^  and 
pya^,  respectively.  Substitution  of  Eqs.  (3)  end  (lb)  into  Eqs.  (10) 
gives  therefore 


*  *py 


l-k  .  2*  do  _  l-k  „  do 

v  *  fpy  #u  *  H  1  “TF5?  ,U  *  dST 

2p/»  *7 


(3-19) 


which  are  the  conponents  of  a  radial  acceleration 

l-k  do 


a  ■ 


n 


Eqs.  (10)  also  furnish  slnple  differential  equations  for  w  and  ▼ 


(3-20) 


PV  *  7  *  7  00t  "  ^T"  oot  *  S’ 


(3*21) 


®v  Z? 1  "  ?  *  V  ®ot  ^  ■ 


l-k  do 
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For  any  vedge  shaped  rag Ion ,  0^  to  0^,  thaaa  equations  daflna  tha 
valocltiaa  in  taraa  of  tbosa  at  ona  aida,  a ay  0^.  Tha  faot  that  tha  acealara- 
tlon  ia  puraiy  radial  suggests  introduction  of  polar  valocltiaa  up  and  u^, 

Ona  than  flnda  tvo  aaparata  equations 

duj 

U  ■  m  I. 


»  d0 


( 3-22) 


(V  sin  0  *  u^) 


1-k  dr 
20  60 


Tha  fact  that  tha  tana  (V  sin  0  ♦  u^)  a ay  rani ah  indlcataa  tvo  distinct 
ranges  of  solutions.  If  tha  valocltiaa  u  and  v  are  saall  •  V"1  and  ▼  «  1, 
Eqs,  (21)  may  be  solved  directlyi 


u(a) 


l-k 

5p7 


f 


cot  0 


da 


60  ♦  u(0o) 


v(0)  -  |=v  l«(0)  -  «(0O)J  ♦  t(0o) 


(3-23) 


Com*  2»  q  ■  0. 

If  the  conditions  (2)  concerning  the  stresses  are  satisfied,  but  e  "  0, 

ju  dY  A 

the  region  la  on  the  verge  of  slip.  In  this  csss,  Eqs.  (11)  give  ^  ‘  '* 

and  tha  right  hand  aides  of  Eqs,  (10)  vanish,  leaving  the  aquations  of 
atatle  equilibrium 


•  0 


dT 

dC 


St 

4? 


•  0 


(3-2l») 
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Substitution  of  Eqs.  (3)  furnishes  tvo  non- linear  simultaneous  differ* 
ential  equations  for  0  and  6 

sin  2d  •  cos  2d)  *  ~  ♦  [f  00s  2d  ♦  sin  2d  )[e(l*k)  ♦  s]  ^  «•  0 

(3-25) 

j±=^ain  2d  ♦  t  coe  2d)  -  (  ♦  jo  os  2d  -  (  sin  2dj  Je(l-k)  ♦  sj  |*  •  0 

These  equations  oan  be  separated,  giving 


Consider  next  the  singular  oase  when  the  determinant  (27)  vanishes. 
Expanding  and  noting  Eq,  (13) ,  one  finds  1 

d  •  0  ♦  T  (3-29) 

•)  Eqs.  (28)  describe  a  region  of  uniform  stress,  the  stresees  satisfying 
the  slip  condition  (2). 
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vhere  the  angle  y  1»  defined  by 


coa  2y 


1-k 

T+3T 


Substitution  of  Eq.  (29)  into  Eq.  (25)  finally  gives 


(3-30) 


do 

d a 


1-k 

TT® 


(3-31) 


ud 


0  •  Ca 


(3-32) 


Except  for  tha  valua  of  tha  exponent,  this  raault  and  Eq.  (18)  art  similar. 
Tha  solution  (32)  aatiafiaa  atatlo  equilibriua  and  is  tha  ms  obtalnsd  in  tha 
classical  atatic  problem  [Raf.  2]. 


B.  Material  with  D1 latency. 

Tha  components  of  tha  strain  rata  tensor,  s  ,  c  and  *  can  bo  sxprossad 

**  » #  *'r 

in  tarns  of  tha  principal  strains,  and  therefore  from  Eq.(2-l8)  by  the  arbitrary 
function  1.  Frocaading  aa  in  Section  A.,  one  finds  that  Kqo.  (10)  ranaln  valid, 
but  in  lieu  of  Eqs.  (11)  one  finds 


*  5?  ’  -  ?  0M  21  *  I?f 

-  d  •  -  ly  sin  28  (3-33) 

Elimination  of  (Ay)  furnishes  the  relation 

(|  -  0  cos  2#  -  2  sin  21  ♦  ♦  j)  -  0 

and  after  substitution  of  Eq.  (13)  one  finds 


•  •  9  ♦  Y 
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( 3—3** ) 


where  y  hu  the  value  defined  in  Eq.  (30).  Proceeding  by  substitution  of 

Eqs.  (3)  into  Eqs.  (10),  one  obtains  a  eet  of  two  eiaultaneoua  non-honogeneous 

J  do  3  d# 

linear  equations  in  ^  and 


sin  26  -  cos  28)  -  ^  ♦  [t  e0«  29  ♦  ,lB  29]  [«U-k)  ♦  ■]  $|  ■  -pV  -  $  ♦  C  7) 

p=Nsin  29  ♦  {  cos  26)  -  C  ~  ♦  [cos  26  -  «  sin  2«][o(l-k)  ♦  »]  $  ■  -pV  -  $  ♦  C  7) 

(3-35) 


Observing  that  the  value  •  given  by  Eq.  (3k)  is  Just  the  one  for  which  the 

deterainant  (27)  of  the  coefficients  of  the  above  equations  vanishes,  on*  must 

conclude  that  these  equations  have  solutions  only  if  a  determinant  formed  by  the 
d6 

coefficients  of  and  by  the  right  hand  sldee  vanishes.  Ignoring  trivial  factors, 
it  is  necessary  that 


(V 


u  ♦  *v) 


C  cos  26  ♦  sin  26 
cos  26  -  (  sin  26 


du 

d{ 


•  0 


dv 

df 


(3-36) 


for  6  ■  0  ♦  y.  The  deterainant  vanishes  if  0  ■  0  or  «  (which  are  trivial  roots) 
or  if  X  *  0.  The  factor  V  -  u  ♦  lr  can  also  vanish,  giving  a  velocity  distribution 
u  ■  V,  v  ■  0,  which  in  tum  again  leads  to  1  •  0.  One  finds  therefore  that  the 
only  possible  notion  oocurs  without  slip.  The  velocities  must  be  constant 


u  ■  constant  v  ■  constant  (3-37) 


while  the  stresses  satisfy  the  equations  of  equillbriun,  and  are  therefore  the 
stresses  for  the  etatlc  problems  treated  in  (2)  and  derived  above  in  Eq.  (32), 


19 


C.  Regions  without  Slip. 


In  locations  without  slip,  whan  tha  inequality  (2-5)  is  satisfied,  tba 
two  materials  considered  may  be  treated  Jointly.  If  no  slip  occurs ,  tbs  vslocities 
u  and  v  in  E qa.  (1,10)  are  necessarily  constants,  such  that  their  derivatives 
vanish,  and  one  is  left  with  the  equation  of  equilibrium 

dox  dt 

dT  '  ‘  4?  *  0 

dT  dOv 

(J.J8) 

In  analogy  with  the  solution  in  [2]  for  the  statio  case,  tbs  construction 
of  solutions  in  the  present  problem  will  involve  the  matching  of  sectors  with 
and  without  slip,  (Fig.  7),  Prior  to  undertaking  the  matching  of  sectors, 
solutions  applicable  to  the  sectors  without  slip  must  therefore  be  obtained.  The 
problem  to  be  studied  concerns  tbs  possible  states  of  equilibrium  in  wedges  of 
arbitrary  angles  d,(Flg.  8a  or  b) ,  loaded  on  the  horizontal  surface  by  a  normal 
load  p(  on  the  inclined  face,  matching  with  the  solution  in  tbs  region  with 
slip  leads  +o  two  conditional 

a. )  The  major  principal  stress  must  have  the  direction  y  required  by  the 
solution  in  A  or  B  above. 

b. )  The  principal  stresses  at  the  boundary  must  satisfy  the  allp  conditions 
Eq.  (2-6). 

However,  there  is  one  additional  condition. 

o.)  Among  the  solutions  of  Eqs.  (38)  only  those  are  acceptable  which  satisfy, 
in  addition,  the  inequality  (2-3),  which  limits  the  possible  state  of  stress  in 
the  Coulomb  material. 
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Eqs,  (38)  can  b«  satisfied  by  expressing  the  stresses  in  the  usual  Mannar 
by  an  Alxy  at rasa-function  P.  In  tha  absaaca  of  a  compatibility  equation, 
thara  is  no  furthar  condition  on  F(0»  and  It  is  laft  arbitrary.  In  thaaa 
circumstances,  families  of  solutions  satisfying  all  requirements  exist,  at  least 
for  certain  ranges  of  0. 

The  multitude  of  solutions  can  be  restricted  by  noting  that  the  assumption 
of  an  elastically  un deformable  material  vas  introduced  because,  in  case  of  slip, 
for  sufficiently  large  values  of  E  and  0,  there  vlll  be  important  situations 
vhere  the  elastic  deformations  can  be  expected  to  be  small  versus  the  slip 
deformations.  However,  in  regions  without  slip,  the  elastic  deformations  ought 
not  to  be  ignored  because  they  are  not  small  versus  the  (non-existent)  slip 
deformations.  Retaining  the  elastic  deformations,  the  acceleration  terms  on  the 
right  hand  side  of  Eqs.  (10)  remain,  and  an  additional  eonditlon  expressing 
elastic  compatibility  applies.  The  appropriate  differential  equations  which 
contain  the  velocity  V  as  parameter,  would  be  those  of  reference  [l],  to  be 
applied  to  a  wedge,  (Fig.  8a,  b) ,  in  lieu  of  the  half  space.  However,  having 
Ignored  the  elastic  terms  required  for  wave  propagation  in  the  slip  region,  the 
analysis  here  is  restricted  to  values 

V  «  d#  «  Op  (3-39) 

For  the  present  purpose,  it  suffices  therefore  to  study  the  solutions  of 
the  elastic  problem  for  the  limiting  ease  V  ♦  0.  In  this  oase,  the  acceleration 
terms  vanish  again,  such  that  the  equations  of  equilibrium  remain  valid|  however, 
the  eompatleility  equation  has  still  to  be  satisfied.  All  stresses  being  functions 
of  «  ■  J,  the  latter  becomes 

(•E  *  9  )  •  0  (3-*»0) 
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The  solution  of  Eqs.  (38  and  bo)  subject  to  tins  boundary  conditions  a.)  and 
b.)  above  is  routine,  but  contains  three  parameters  0,  k,  and  s.  Having  obtained 
the  solution,  it  remains  to  be  determined  for  what  range  of  the  parameters 
condition  0.)  above  is  satisfied.  This  complicates  the  presentation  appreciably 
and  the  entire  matter  is  therefore  relegated  to  Appendix  A{  only  pertinent 
conclusions  are  listed  belovt 

1,  For  the  "incompressible"  material,  when  at  the  inclined  face  of  the  vedge 
Y  ■  ±  j“,  it  vas  found  that  solutions  satisfying  all  conditions  (a,  b,  and  0) 
exist  in  the  range 

fiBiys  (S-M) 

provided  the  directions  of  the  major  principal  stress  is  as  sbovn  in  (Fig.  9a  or  b), 
1,  e.,  when  <fx  makes  the  angle  ^  ♦  0  with  the  horizontal  side  of  the  wedge.  For 
this  orientation  of  the  major  principal  stress,  no  solution  exists  for  other 
angles.  Just  one  solution  exists  for  each  combination  of  k  and  s.  The  value  of  • 
on  the  Inclined  face  of  the  wedge  is  given  in  Eq.  ( A<-35 )  •  For  the  wedge  angles 
S  ■  J  and  the  solutions  represent  uniform  states  of  stress,  and  are  of  the 
type  previously  found,  Eq,  (28),  for  the  special  ease  e  ■  0. 

2.  For  the  "incompressible"  material,  similar  solutions  exist  if  the  direction 
of  the  major  principal  stress  with  the  horicontal  is  (I  •  ^  ),  (Fig.  10a,  b),but 
only  if 

B  ■  ^  or  6  •  (3-*»2) 

For  these  two  values  of  B,  a  solution  exists  for  any  combination  of 
k  and  s.  These  solutions  represent  states  of  uniform  stress,  and  are  of  the 
type  previously  found  in  Eq.  (28), 
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3.  In  part  A  of  this  Motion,  under  Cm*  2,  a  aolution  without  slip  with 
a  non-uniform  state  of  stress,  Eqs.  (29  to  32),  has  been  developed  for  the 
"incompressible"  material.  In  view  of  our  present  reasoning,  the  compatibility 
condition  (ho)  should  be  satisfied  in  a  region  without  slip.  The  solution 
Eqs,  (29  to  32)  must,  therefore,  be  checked  in  this  respect.  Rot  unexpectedly, 
one  finds  that  compatibility  is  violated,  such  that  this  solution  will  not  be 
considered  further. 

h.  In  the  case  of  the  material  with  41 latency  when  the  angle  y(l)  satisfies 
the  equation  cos  2y  ■  solutions  also  exist  for  ranges  of  vedge  angles  f . 

The  details  are  not  presented,  for  reasons  discussed  in  Section  3  where  the 
material  with  dilatancy  is  discussed. 

3.  For  completeness  sake,  the  possibility  of  regions  of  no  slip  between  two 
regions  with  slip  has  been  considered,  beoause  this  possibility  would  permit 
geometries  different  from  (Fig,  7).  Such  regions  exist,  provided  |  •  j.  The 
solutiona  are  uniform  states  of  stress,  covered  by  Eq.  (28).  (nte  situation, 
however,  oan  not  be  utilised  in  the  construction  of  solutions.) 

6.  It  is  of  Interest  to  know  for  what  values  of  p,  end  up  to  what  values  of 
p,  elastic  solutions  for  the  half  space  exist  without  slip,  1.  e.,  solutions  where 
the  inequality  (2-3)  is  satisfied  everywhere.  Eqs.  (38,  ho)  apply  again  with 
appropriate  boundary  conditions  at  0  •  0,  i,  (Fig.  12),  and  are  considered  la 
Appendix  A. 

If  there  is  no  cohesion,  s  •  0,  no  solution  satisfying  Iq.  (2-3)  exists 
for  p  4  0. 

However,  if  s  4  0,  an  unusual  situation  occurs.  Ths  solutions  of  the 
differential  equations  for  the  stresses  satisfying  the  boundary  condition  on  the 
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surface  art  not  unique,  as  a  stats  of  uniform  horisontal  stress  9(  ■  constant 
can  be  added  to  any  solution  one  may  find.  Ibis  arbitrary  horisontal  stress 
enters  the  condition  for  slip,  Eq.  (2-3),  and  the  value  of  p  up  to  which  no 
slip  occurs  depends  on  the  assumed  intensity  of  the  horizontal  state  of  stress. 

There  is,  however,  a  most  "beneficial"  state  for  which  the  load  p  becomes  a 
maximum.  This  "distinguished  solution"  defines  the  value  of  p  above  which  solutions 
without  slip  can  not  possibly  exist,  A  plot  of  these  critical  values  of  p  vs,  k 
is  given  in  Appendix  At(fig.  A. 4). 
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k.  SOLUTIONS  FOR  THE  "INCOMPRESSIBLE"  MATERIAL. 
a.  Construction  of  Solutions. 

In  this  Motion,  the  solutions  applying  for  ths  "inconpresslble”  material. 
Csss  A  in  Ssotlon  2,  will  bs  obtained.  lbs  svsilsbls  lngrsdlsnts  srsi 

1.  A  solution  in  s  rsglon  with  slip,  trsatsd  undsr  A,  Csss  1,  la  tbs 
prsrious  ssotlon,  It  vss  found  tbst  la  sucb  a  rsglon,  tbs  najor  principal 
stress  aakss  an  angle  ♦  ^  with  tbs  position  rector  0,  If.  (3-1*0,  while  s  Is  glrsa 
by  Iq.  (3-18).  It  will  be  seen  that  tbs  ease  of  tbs  lower  sign,  (Fig,  13),  will 
bs  used.  In  this  ease 

s  •  Ce*  2  TST  0  .  jJ-  (*•-!> 

whs re  C  Is  an  open  but  positive  constant.  Tbs  equation  Indicates  that  S  Increases 
wltb  Increasing  values  of  0, 

2.  Of  solutions  without  slip  obtained  under  A,  Case  2  In  Section  3,  only 
the  ones  defined  by  Kq,  (3-28) 

9  •  oonst,  •  •  const,  (4-2) 

representing  a  state  of  uni font  stress  oan  be  considered.  The  other  solution, 

Kqs.  (3-29  to  32)  is  not  to  be  used,  for  reasons  stated  in  Section  3  under  C, 
iten  3,  Solutions  of  tbs  type  of  Iq,  (2)  will  be  used  for  45°  wedges  as  shown  In 
(Fig.  14a,  b).  In  (Fig,  l4a),  tbe  najor  principal  stress  is  Inclined  at  g  ■  j  , 
while  In  (Fig.  l4b)  tbe  inclination  Is  •  ■  0, 

3.  Additional  solutions  without  slip  is  wedge  shaped  regions  are  described 
In  Section  3  under  C,  lten  1.  Of  Interest  are  solutions  for  wedges  |i  |  < 

(Fig.  19). 
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Let  u>  bov  consider  the  problem  of  o  stop  load  progressing  with  velocity 


V,  (rig.  2),  for  Mtorlal  proportloa  k  aad  at  for  various  valuta  of  tha  load) 
starting  vlth  aaall  valuta  of  For  sufficiently  snail  valuta  of  tbia  ratio. 

It  bat  bats  shown  la  Appendix  A  that  (tlaatlo)  solutions  without  slip  salat, 
provldtd  that  p  la  lass  thaa  a  critical  value  p(,  given  In  (Fig.  A-U).  As 
pointed  out  In  Section  3  under  C,  Itca  6,  the  solutions  are  not  even  unique, 

1.  e«,  there  are  fanllles  of  solutions  If  p  <  p^,  but  there  la  Just  one  If  p  *  p#. 
The  erltloal  value  p^  vanishes  If  a  •  0.  If  cohesion  la  present,  assuming  likely 
values  of  k,  the  Unit  p§  is  2a  to  3s,  1.  a.,  at  an  uninterestingly  lot/  level  of 
stress. 


When  constructing  solutions  vlth  slip,  It  la  clear  that  solution  1.  above 
Is  Insufficient  to  satisfy  the  condition  at  either  the  loaded  or  the  free  surfaoe, 
beoause  the  directions  of  the  principal  stresses  at  the  surface  eannot  be  at 
angles  ♦  Jto  a  surface.  Transition  regions  without  slip  are  necessary  as  Indicated 
In  (Fig.  7). 

The  solutions  1,  2  end  3  shove  can  be  ooa&lned,  (Fig,  16),  to  fora  a  eon* 
tlnuous  solution  for  the  half  space,  consisting  of  a  non-slip  wedge  of  angle 

V  If  f 

^  <  I  <  on  the  loaded  side,  a  non-slip  wedge  of  angle  f  on  the  unloaded  side, 
and  an  Interior  wedge  of  angle  (^  -  |)  in  which  slip  occurs. 

Consider  first  the  cue  I  ■  such  that  slip  oocurs  In  n  90°  wedge.  Under 
the  load,  the  state  of  stress,  (Fig,  lHa),  Is  uniform,  the  aajor  principal  stress 
being  vertical!  In  this  range 

«  ■  p  •  e  (*-3) 

In  the  adjoining  sector  with  slip,  the  major  principal  stress  nakes  therefore 


26 


the  angle  -  ^  with  the  poeltlon  rector \  obtaining  the  value  of  C  in  Eq.  (1)  by 
■etching  the  values  of  •  at  j  ■  one  finds 

J.d.JI,  ,  .  (p  ♦  g_).-  2  &  •  «>.  (WO 

As  the  direction  of  the  major  principal  stress  in  the  region  0  «  ^  aust 
■etch  the  direction  in  the  slip  region,  o^  Bust  be  horlsontal,  such  that  (Fig,  ihb) 
applies,  and  one  finds  that  thd  solution  described  requires  a  vertical  preeeure 

at  0  •  Oi 


P 


o 


(M) 


The  surface  at  0  •  0  is,  however,  a  free  surface,  where  the  pressure  should 

vanish,  p  •  0.  Solving  Eq.  (5)  for  the  value  p  ip  for  vhleb  p  ■  0,  one 
O  L  v 

obtains 

"i  ‘  ☆[.  ra  * .  k]  <l"<> 

such  that  a  steady  state  solution  for  the  speelfie  value  p  ■  p^  has  been  obtained* 

The  subeoript  L  has  been  selected  because  this  value  represents  a  Uniting  situation, 
as  will  be  recognised  later* 

One  oould  now  proceed  to  vary  the  angle  |  in  (Fig,  1 6)  sad  try  to  obtain 

solutions  for  values  p  0  p  .  However,  the  presence  of  cohesion  masks  the  true 

L 

state  of  affairs,  and  as  a  first  step  only  the  case  a  •  0  will  be  considered.  In 

this  ease,  Eqs.  (1,  k,  3)  retain  the  exponential  terns  on  the  right  hand  side  only, 

and  lq.  (6)  gives  p  ■  0,  1,  s,,  only  a  trivial  solution,  for  p  •  0,  has  been 
L 

dets  mined. 


If  one  uses  different  angles  I  in  (Fig.  16)  on  the  loaded  side,  one  agala 
finds  only  trivial  solutions  p  ■  0,  as  is  easily  confirmed.  The  relation  between 
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p  sad  o  on  the  inclined  face  for  the  wedge,  (fig.  15).  is 


m  m  2p  (1— cos  26)  -  s( 1-cos  26  ♦  tin  28  -  28  coo  2fl) 

Uek)  (1-cos  20)  ♦  (1-k)  (.in  20  -  28  coo  20)  (M) 

When  •  •  0(  the  value  of  a  at  tha  inclined  face  is  therefore  proportional 
to  pi  aatehing  of  the  solution  (1)  gives  C,  and  the  stress  e  at  0  ■  both 
proportional  to  p.  This  is  to  be  matched  with  the  l»5°  wedge  shown  in  (Fig,  lkb). 

If  s  ■  0(  and  the  stress  at  0  ■  0  vanishes,  the  stress  o  in  this  wedge  necessarily 
vanishes ,  and  p  Bust  also  vanish. 

There  is,  therefore,  for  s  ■  0,  no  solution  for  the  boundary  value  problem 
formulated,  i.  e.,  for  a  finite  value  of  p  behind  the  moving  front,  and  no 
pressure  ahead  of  the  front.  From  a  purely  mathematical  point  of  view,  this  is 
connected  with  the  fact  that  the  equations  are  hyperbolic,  and  saving  prescribed 
conditions  on  one  side,  one  cannot  simultaneously  prescribe  conditions  on  the 
other.  However,  the  physical  problem  must  have  a  solution  -  there  must  be  a 
response  -  if  a  pressure  wave  is  applied  to  the  half-space.  The  mathematical 
and  physical  concepts  can  be  reconciled,  by  allowing  for  the  possibility  of 
particles  of  the  material  being  expelled  from  the  unloaded  surface  ahead  of  the 
pressure  front.  It  is  reasoned,  that  the  non-existence  of  eolutlone  of  the  problem 
for  certain  applied  loads  in  the  formulation  used  so  far  indicates  that  the  body 
cannot  exist  under  these  loads,  and  will  disintegrate.  The  process  of  disintegra¬ 
tion  has  been  studied  for  a  simple  equivalent  situation  in  Appendix  B  on  a  mechanical 
model  exhibiting  Coulomb  properties.  A  mechanism  is  demonstrated  which  expels > 
partlclee  on  an  unloaded  surface,  the  momentum  of  the  particles  providing  a 
reactive  prceeure  required  to  obtain  a  solution  in  the  Interior  of  the  body. 

In  thle  revieed  situation,  regardless  of  the  value  of  s,  Kqs.  (|,  k)  represent 
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•  solution  of  the  steady  state  problem  for  p  _>  p  ,  but  a  reactive  pressure 
according  to  Eq.  (5)  is  exerted  by  particles  expelled  from  the  surface  ahead  of 
the  pressure  front  p. 

Continuing  the  discussion  for  s  ■  0,  the  limit  p  ■  0,  and  the  solution 

L 

Eqs,  (It,  5)  applies  for  all  values  of  p.  However,  it  is  not  the  only  solution. 
Varying  the  angle  I  of  the  non-slip  region  in  (Fig.  16),  other  solutions  can  be 
obtained.  The  value  of  e  at  0  ■  «  -  0  follows  from  Eq.  (7)  for  s  ■  0: 


o(*  “  "  \T+k) . ri-coi  2^  MMtrfii'nTi  7  2*  cos  2*) 


(4-8) 


In  the  adjoining  slip  region  o  is  given  by  Eq.  (l)fvhere  C  is  obtained  by 
matching  values  at  0  ■  t  -  I.  In  the  non-slip  region,  0  <  ((  <  j,  the  value  of  • 
is  constant,  and  one  finds 


0(0)  ■  o(«  -  0)  e"2  1+k  *  •) 


(4-9) 


The  required  reactive  pressure  to  be  exerted  by  expelled  particles  at 
0  ■  0  becomes 

1-k  /3t  \ 

Pq  •  ke(0)  •  ko(w  -  I)  •  T*k  UT  *  *)  (4-10) 

In  addition  to  the  solution  for  I  ■  45°  found  first,  there  exists  therefors 
a  whole  family  of  solutions  for  J  <  I  <  All  these  solutions  are  mathematically 
unobjectionable |  this  lack  of  uniqueness  in  steady-state  problems  was  discussed 
in  the  Introduction  and  is  nothing  unusual. 

Comparing  the  possible  solutions  for  various  angles  (Tig,  16),  it  is 
seen  that  the  one  for  |  ■  45°  where  slip  ocours  in  a  90®  wedge  has  distinctive 
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properties,  end  it  will  be  called  the  Nprincipal  solution".  Among  all  solutions 

for  a  given  value  of  the  applied  pressure  p,  it  requires  the  smallest  pressure  p^ 

exerted  by  the  particles  expelled  on  the  free  surface.  To  demonstrate  this 

consider  (Fig.  17).  It  above  the  values  of  a  for  the  "principal  solution"  as 

v  3* 

function  of  0  plotted  as  Curve  It  in  the  range  £  <  0  <  the  curve  decays 
exponentially  as  given  by  Eq.  (It),  while  the  remainder  is  horizontal.  Curve  II 
in  (Fig.  17)  is  the  locus  of  the  values  o(t  .  6)  according  to  Eq.  (8),  i.  e . , 

the  value  of  o  at  the  Interface  of  the  slip  and  non-alip  regions.  The  and  points 

1  w 

of  the  curve  are  unity  and  J  respectively,  and  it  has  a  minimum  at  B  ■  as  can 

be  found  by  differentiation  of  Eq.  (8).  One  can  easily  demonstrate  that  between 

w  3*  «  v 

■JJ"  and  ]J“  Curve  II  is  always  above  Curve  I,  and  between  |^  and  jjT  the  slope  of 

Curve  I  la  always  steeper  (downward)  than  that  of  Curve  II. 

Curve  III  indicates  the  situation  for  an  intermediate  value  of  |.  At  the 
Interface  0  ■  («  -  ft)  the  value  of  o  is  defined  by  point  C(  the  value  a  then 
decreases  exponentially,  1.  e.,  the  decay  curve  C-D  is  a  portion  C*»D*  of  the 
curve  A  B  displaced  horizontally,  euch  that  the  final  value  at  the  surface, 
q(0),  and  pQ  a  ko(0)  is  necessarily  larger  than  for  the  "principal  solution". 

In  the  general  case,  vhen  cohesion  does  not  vanish,  the  situation  is  more 

complicated.  For  low  pressures,  p  <  p#l  there  are  elastic  solutions  without 

slip  anywherej  for  p  >  p  there  la  tha  "principal  solution",  requiring  disintegra* 

L 

tlon  at  the  free  surface  if  p  >  p  ,  while  for  p  ■  p  no  disintegration  occurs. 

L  « 

(Fig.  18)  shows  a  plot  of  p  and  p.  versus  k,  indicating  that  everywhere  p.  *  p_. 

I  L  L  I 

There  is,  therefore,  a  gap  between  tha  levels  of  p  and  p.  where  neither  the  elastic 
nor  the  "principal  solution"  applies.  For  such  values  of  p  solutions  can  be  eon* 
structed  by  utilizing  wedge  angles  I  »  (Fig.  16).  If  s  »*  0  the  Curve  II,  (Fig.  lT)t 
describing  a  at  the  inclined  faoe  is  modified,  but  the  general  situation  remains 
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similar,  such  that  for  any  value  p  >  p  >  p  ,  a  value  of  $  exists  for  vhlob  p 

L  s  0 

given  by  Eq.  (10)  will  vanish.  In  other  words,  In  the  range 

p  >  p  >  p  ( U-ll) 

L  • 

the  Eqs.  (8,  9,  and  10)  and  the  condition  PQ  ■  0  define  a  "principal  solution" 
in  which  slip  occurs  in  the  range  £  «  0  <  «  -  0,  where  j[  <  t  <  (Fig,  16). 

For  this  solution  the  pressure  on  the  surface  ahead  of  the  wave  vanishes  such 
that  no  disintegration  occurs.  (It  would  also  be  possible  to  find  polutlons  in 
this  range  where  pQ  >  0,  requiring  disintegration.  Such  solutions  are  not 
considered, because  a  "principal  solution"  with  PQ  ■  0  exists.) 

The  Method  of  selecting  the  "principal  solution"  is  not  very  satisfactory 
as  intuitive  reasoning  has  to  be  employed.  More  definite  conclusions  showing 
that  these  solutions  are  the  ones  approached  in  transient  cases  after  a  long  tine, 
would  require  consideration  of  the  starting  conditions ,  which  would  pose  a  major 
yet  unsolved  problem.  As  an  alternative,  one  might  study  the  stability  of  the 
various  solutions  against  perturbations.  This  might  permit  elimination  of  some 
excess  solutions,  but  such  an  approach  might  not  be  successful.  Stability 
considerations  are  unable  to  resolve  the  ln«k  of  uniqueness  in  the  elastic  rolutloms 

for  pep. 

s 

An  additional  point  is  worth  discussing  and  illustrates  the  inherent 
complexities  of  the  type  of  problem  considered.  (Fig.  18)  shows  plots  of  p  sad  p. , 

S  L 

and  also  of  the  value  of  a  statio  load  p  which  can  be  supported  by  a  half  space 

P 

according  to  the  classical  analysis  by  Prandtl  [2],  It  is  seen  that  p  »  p  , 

r  L 

such  that  the  present  paper  contains  in  the  limit  V  ♦  0,  solutions  where 
steady-state  motion  occurs  at  pressure  levels  below  p^  (where  equilibrium  is 
possible).  It  is  outside  the  province  of  this  paper  to  investigate  if  this  in* 
dleates  an  instability  of  the  static  solution  (2)  for  the  material  considered 
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here  or  not.  It  in  important ,  however,  to  point  out  that  tba  fact  that  tha 
limit  V  ♦  0  of  tha  praaant  aolution  doaa  not  give  tha  atatie  altuatlon  ia  not 
dua  to  a  aathamatieal  arror.  In  tha  atatie  oaae,  vban  tha  load  p  »  p?  la 
applied,  allp  muat  occur  in  boom  locationa  bafora  tha  atata  of  equilibrium 
ia  raacbed.  If  tha  load  la  moved  with  a  velocity  V,  however .email,  alip  haa 
to  occur  ooatinuoualy,  and  a  aliplaaa  aolution  at  tha  level  p  •  p_  la  therefore 
not  poaalble.  Tha  dynamie  aolution  for  V  ♦  0  doaa  not  approach  tha  atatie 
aolution. 
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b«  Discussion  of  Strains  &.b4  Accelerations.  due  to  the  "Principal  Solutloq" 

In  the  preceding  subsection  (a)  solution*  with  slip  burs  bssn  obtained 
covering  ths  rang*  of  prsssur«st  p  >pj(  vhsrs  p#  dsflnss  tbs  Unit  up  to  which 
elastic  solutions  without  slip  exist.  Xf  oohesion  is  present,  s  4  0,  then  there 
is  a  narrow  range  p,  <  p  t,  where  the  "principal  solution  "  selected  involves 
slip  hut  does  not  require  surface  disintegration  ahead  of  the  pressure  pulse.  For 
larger  pressures,  p  >  p^,  solutions  exist  only  with  slip  and  surface  disintegrate 
tioo,  the  "principal  solution"  being  the  one  with  the  least  reactive  pressure  due 
to  expelled  particles. 

It  is  naturally  of  Interest  to  discuss  the  magnitude  of  accelerations  and 

strains  associated  with  the  principal  solutions.  For  these  quantities ,  there 

being  only  small,  quantitative  differences  between  the  reaves  p  <  p.  and  p  >  p  , 

I*  L 

only  the  latter  will  be  considered. 


A  quantity  of  interest  is  the  total  permanent  strain  produced  at  any  point 
due  to  the  passing  of  the  shock  wave  on  the  surface.  This  strain  gives  an  indica¬ 
tion  of  the  deformation  to  whieh  a  target,  say  a  cylindrical  shell.  (Fig,  19). 
would  be  subjected  if  the  presence  of  the  target  would  not  affect  the  force  field*). 

To  obtain  the  total  strain,  the  strain  rates  at  the  target  at  depth  y  must  be 
obtained,  and  Integrated  with  respect  to  time.  Substituting  Xq.  (>-23)  into  Kq.  (3-11) 
and  noting  that  the  angle  t  in  the  solution  employed  is  defined  by  Sq.  (3-lU)  with 
the  lower  sign,  one  finds 


e  « 


1-k  d«J 

Wy  4? 


(>*-12) 


This  expression  applies  in  the  slip  region  ^  <  0  <  jp’.  The  value  of  •  is 

*T TRwfl&er  or  not  the  target  affects  the  free  field  substantially  or  not  depends 
on  relative  stiffness,  but  little  is  known  about  the  natter  at  present* 
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given  by  Eq,  (1*) , 


(4-13) 


which  pemite  the  determination  of  the  strain  rates  and  from 

Eqs.  (3-6). 

To  integrate  the  resulting  expression  with  respect  to  time,  it  is  noted 

that 


dt  5 


V  sin20 


d i 


Hence 


•  v  jv  “ 

•w  '  |«: 


2A 


dt  ■  -  A 


3w 

■r 


f  t « „ 

,  .zw"  cot  0  Slf 

IT 


«  (1  -  cot2*)  e2 


where 


(«*-H0 


(**-15) 


(i-16) 


A  ■  isSw  (n  iliL  45iL.\.'2  (!♦*) 

2^V?  V  T3T  4  14k/* 

e  .  3s 

The  limits  in  the  above  integrals  are  those  of  the  slip  region  ^  <  9  <  £"% 
The  nature  of  the  Integrals  requires  numerical  integration.  Another  quantity  of 
Interest  is  the  acceleration  history.  The  accelerations  in  the  non-slip  regions* 
(Fig,  19),  vanish)  in  the  slip  region,  ^  (3-20)  defines  the  radial 

acceleration.  Using  Eq,  (13) 


a  ■ 


; *  8s)  -2cot'1(?f 


(«wl7) 


Example i  The  phyaieal  meaning  of  these  results  is  best  discussed  for  a 
typical  case.  We  select  k  ■  ^  and  no  cohesion,  s  ■  0.  Using  Iqs.  (19)  to  compute 
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the  permanent  strain*,  one  finds  the  principal  permanent  etralne 


*1.2  li  "-in  ^7 


(•> 


the  direction  being  defined  by  e  v  0.5.  To  interpret  the  result  it  will  be 
coopered  with  the  magnitude  of  elastic  strains  which  are  of  the  order  £  ,  Hot in* 
the  expression  for  the  velocity  of  P*waves  in  an  elastic  medium  d2  ■ 


p  pTI+w  TOT5T  • 


Eq,  (a)  may  be  written 


c 

1.2 


(b) 


Due  to  the  basic  assumptions  on  which  this  paper  is  based,  V  must  be  small 
versus  c^(  such  that  the  term  in  parentheses  is  much  larger  than  unity.  The 
permanent  strain  is,  therefore,  much  larger  than  elastic  strains  at  this  stress 
level. 


Eq.  (b)  also  permits  an  estimate  of  the  range  of  applicability  of  the 

present  theory,  which  ignores  elastic  strains!  It  should  be  applicable  as  long 

as  the  slip  strains,  Eq.  (b),  are  appreciably  larger  than  elastic  ones,  This 

V 

condition  is  reasonably  satisfied  up  to  Mach  numbers  rr-  <  0.2,  and  the  theory 

**  y 

presented  will  therefore  be  applicable  up  to  such  values  of  . 

*P 

The  value  of  the  radial  acceleration  a,  given  by  Eq.  (17)  is  plotted  in 
(Fig.  20).  The  acceleration  a  vanishes  if  |y|  >  1|  its  maxima  value  in  the 
range  J^|  «  1  is 

max  a  \  0.24 

-  py  vo; 

For  comparison  purposes  one  can  express  p  in  terms  of  the  velocity  of 
elastic  shear  waves  «  end  of  the  modulus  of  rigidity,  0, 

^2 

max  a  \  0.24  (d) 
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This  value  may  be  compared  with  the  acceleration  due  to  a  step  vave  In 
the  elastic  case,  which  can  he  obtained  from  the  results  In  Ref.  [l].  Using 
the  curve  in  Fig.  5  of  this  paper,  the  effect  of  a  step  vave  can  be  found  by 


integration.  It  is  seen  that  the  acceleration  given  by  &<i.  (d)  is  of  equal 

V 

magnitude  as  the  elastic  one  for  M  ■  0.5  which  corresponds  to  -—  ■  0.25.  Aa 

T  °p 

the  accelerations  in  the  elastic  case  decrease  with  the  fourth  power  of  V,  the 

present  solution  gives  again  substantially  larger  values  than  [Ref.  1]  provided 

~  <  0.2. 

c 

P 


It  is  important  to  note  that  the  accelerations  do  not  depend  on  the  velocity 

>2 

V,  but  that  the  total  strains  are  proportional  to  V  .  The  Independence  of  a  is 

due  to  the  fact  that  the  accelerations  are  solely  determined  by  the  streee  field, 

by  virtue  of  Eqs.  (3-1),  and  the  velocity  V  does  not  explicitly  appear  in  these 

equations.  To  find  strain  rates  and  strains,  one  and  two  integrations,  respectively, 

with  respect  to  time  are  required.  Because  of  Eq.  (1*0,  each  such  Integration 
1 

Introduces  a  factor  7,  explaining  the  fact  that  the  strains  are  proportional  to 


The  above  comparison  of  effects  permits  an  interesting  prediction  for  the 

V 

range  of  larger  Mach  numbers,  ~  >  0.2,  for  a  material  having  combined  elastie  and 

P 

Coulomb-slip  properties.  While  for  low  Mach  numbers,  the  slip  effects  are  much 

V 

larger  than  the  elastic  ones,  this  is  not  the  case  for  larger  values  of  “.  For 

P 

such  values,  the  occurrence  of  slip  will  not  change  the  order  of  magnitude  of 
acceleration  and  strains  found  by  a  purely  elastic  analysis  ignoring  slip.  One 
can  understand  this  conclusion  by  noting  that  for  large  values  of  V,  the  slip  stage 
lasts  only  a  short  time  such  that  displacements  due  to  slip  have  only  little  time 
to  develop,  and  will  remain  of  the  same  magnitude  as  the  elastic  ones. 
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Another  prediction  for  result!  at  larger  Mach  numbers  concerns  the  dis¬ 
integration  of  the  surface  ahead  of  the  pressure  wave*  This  effect  vlll  occur 
also  at  larger  Mach  numbers,  but  of  course  only  for  sub-selsalc  velocities. 


3 1 


5.  THE  MATERIAL  WITH  DILATANCY 


For  this  mat* rial  a  solution  has  bsen  obtained  in  Section  3  under  B( 
whars  tba  diraotion  of  tha  major  principal  stress  mskaa  sn  angle  +  y  with  the 
radial ,  (Fig,  21) twhara 

cos  2y  •  (5-1) 

This  solution  could  be  matched  with  regions  of  uniform  stress  if  ths 
angles  are  selected  as  shown  in  (Fig.  22).  One  can  then  find  a  value  pi  p^ 
where  the  surface  ahead  is  a  free  surface ,  the  solution  agreeing  with  the  static 
one  found  by  Prandtl  [2],  For  larger  values  of  the  load,  p  >  pp  (  one  obtains 
a  solution  requiring  a  pressure  pQ  ahead  of  the  pressure  front,  pQ  being  supplied 
again  by  the  disintegration  of  the  surface. 

One  could  further  investigate  if  solutions  exist  for  wedges  without  slip 
at  angles  other  than  90°  •y  andy,  utlliaed  in  (Fig,  22).  This  investigation 
has  not  been  mads  for  reasons  explained  hereafter* 

In  8eotlon  3  under  C,  it  vas  concluded  that  ia  regions  without  slip  the 
equations  of  elustle  compatibility  should  be  satisfied,  such  that  the  solutions 
found  in  this  paper  would  apply  for  an  aotual  elastic- Coulomb  material  in  the  limit 
when  E  and  0  are  large.  The  fact  that  elastic  deformation  in  areas  where  slip 
ooours  are  neglected  la  the  present  paper  is  defensible  as  long  as  ths  slip  deforma¬ 
tions  are  large  versus  the  elastio  ones*  This  argument,  unfortunately,  does  not 
apply  for  the  solutions  found  in  Section  3  under  B,  because  it  vas  found  that  ths 
solution  requires  AjbO,  and  according  to  Sq*  (3-33)  all  strains  and  dsformatlons 
vanish*  In  these  circumstances  ignoring  of  the  clastic  dsformatlons  is  not 


3Q 


Justified.*)  if  ou  desires,  ss  ve  do,  to  obtain  s  solution  vhiob  applies  to  a 
material  having  elastle-Coulonb  properties* 

In  these  circumstances,  it  vas  decided  not  to  pursue  all  the  details 
of  the  possible  solutions,  as  vas  done  in  Section  k  for  the  alternative  material* 

It  would  be  unreasonable  to  waste  effort  on  a  solution  because  all  deformations 
vanish  and  one  obtains  in  the  end  rather  trivial  results i  All  accelerations  and 
the  free  field  strains  (corresponding  to  Eq.  (*ul5)  for  the  other  material)  vanish, 
the  only  non«vanlshing  results  concern  stresses. 

For  materials  with  Di latency  an  analysis  Including  elaatlo  deformations 
is  neoessary,  and  is  planned  for  the  future.  Incidentally,  one  can  easily 
understand  that  the  coupling  of  slip  and  volume  expansion  in  Xqs.  (2*9)  requires 
that  a  meaningful  analysis  must  retain  the  effects  of  elastic  volume  changes* 


•)  In'  order  to  avoid  misunderstandings,  it  should  be  stressed  that  the  above 
objection  to  Ignoring  elastic  effects  only  applies  in  the  dynamic  prbblsm 
studied  here.  The  objection  does  not  apply  in  a  static  ease,  treated  by 
Prandtl  (2),  and  the  above  reasoning  should  not  be  misconstrued  into  an 
Objection  to  the  latter  solution. 
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6.  SUMMARY  OF  CONCLUSIONS 


The  possibility  of  ateaty-stat*  aolutiona  with  alip  haa  been  lnvaattgatad 
vbaa  tba  half-space  la  loaded  by  a  moving  atap  load,  (Fig.  2),  Raatricting  the 
atudy  to  lov  Macb  numbers,  tha  alaatic  da format Iona  vara  neglected  In  the  alip 
region*.  Two  matarlala  vara  considered^  both  follov  tha  Coulomb  rule,  defining 
tha  atraaaaa  at  alip,  but  have  different  flow  rule*. 


One  material  la  aaauaad  to  remain  conatant  in  volume  during  alip,  while 

the  other,  ualng  the  concepta  of  Drucker  and  Prager  [2]  exhibits  Di latency.  For 

the  latter  material  only  unaatla factory  redulta  are  obtained,  becauae  alaatic 

deformation*  were  neglected.  For  tha  other  material,  the  reaulta  obtained  are 

V 

expected  to  be  valid  up  to  Mach  number*  -*"•  ~0.2. 


A  major  finding  of  the  preaent  atudy  la  the  recognition  that  for  preaauree 

p  above  a  critical  value,  dlalntegratlon  of  the  aurfaoe  ahead  of  tha  preaaur*  wave 

muat  occur,  producing  a  preouraor  of  duat  particle*.  The  mechaniem  of  the  dia- 

integratlon  haa  been  atudled  on  a  mechanical  modal,  and  it  haa  been  ahown  that 

element*  in  a  boundary  layer  near  the  free  aurfaoe  will  be  accelerated  upward  and 

finally  expelled, (Fig,  23) »  Thla  manner  of  dlalntegratlon  oooure  for  both  type* 

V 

of  material*  oonaidered,  and  ia  expected  to  occur  beyond  the  range  <  0.2  of  the 
preaent  atudy, but  only  for  aubaelamic  velocitiea  V. 


Detailed  reaulta  for  the  material  without  rolvaae  change  during  alip  are 
preaented  in  8eotion  4,  obtaining  atreaaea,  acceleration*  and  permanent  atralna. 
It  la  demonatratad  (Section  4,b)  that  aooelerationa  and  permanent  stralna,  ia  tha 

y 

range  4-- *5*  0.2,  are  aubatantially  larger  than  the  correaponding  quantities  in  an 
P 

elastic  medium  vhen  no  alip  occur*. 
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The  results  can  be  utilised  to  study  the  equivalent  p rob lee  for  e  looking 
nedlua  having  Coulomb  properties  behind  the  looking  front*  Such  on  Investigation 
Is  under  vay. 

For  the  notarial  with  Dilstency  only  stresses  can  be  obtained  as  the 
analysis  degenerates  and  accelerations  and  strains  are  found  to  vanish*  To  obtain 
meaningful  results  elaatle  effects  aunt  be  retained. 
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APPENDIX  A 


STATIC,  ELASTIC  SOLUTIONS  FOR  WEDGE  SHAPED  REGIONS  AND  FOR  THE  HALF-SPACE 
WITHOUT  SLIP. 


As  background  for  Section  3C,  the  elastic  stresses  are  obtained  for 
vedge  shaped  regions, (Pig.  A-l),  Using  for  convenience  polar  coordinates 
and  stress  components  and  i  ,  the  following  conditions  on  the 

surfaces  are  to  be  satisfied. 


On  the  surface  ♦  »  0: 

o#(0)  -  -  P,  t  -  0  ,  (A-l) 

while  on  the  surface  ♦  -  0,  the  direction  of  the  major  principal  stress  shall 
be  in  a  specified  direction,  defined  by  the  angle  y(0),  (Fig.  A-l).  (Of  interest 
will  be  the  angles  7(0)  •  +  ^).  Further,  at  ♦  -  0  the  two  principal  stresses 
e^  and  a2  are  required  to  satisfy  the  slip  condition,  such  that 


a 

e 


1 

2 


(•  ♦  •) 

ks  . 


(A-2) 


Having  obtained  solutions for  these  boundary  conditions,  it  is  to  be  investigated 
which  of  the  solutions,  if  any,  satisfy  the  Coulomb  rule 


-  o2  >  k*  ■  -  k^  ♦  s) 

everywhere . 


(A-3) 


The  equations  of  equilibrium ,( 3- 30  \  in  polar  coordinates  for  stresses  which 
are  Independent  of  r,  become 


+  Cj,  -  ■  0 

d 

ar +  2t  • 0 


(A-4) 
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while  the  compatibility  equation  (3-l»0)  become* 

d2 

— 5  +  ®  (A-5) 

d  ♦ 

The  general  solution  of  the  differential  equations  (4,$)  may  b*  written 

s  -  A  +  B4  -  C  sin  24  -  D  cos  2# 
r 

■  A  ♦  B4  +  C  sin  2#  +  D  cos  24  (A-6) 

T  ■  -  5  .  (j  cot  24  +  D  sin  24 

where  A,  B,  C  and  D  are  arbitrary  constants.  Determining  A  and  B  from  the 
boundary  conditions  (1)  gives  the  stresses  in  terms  of  C  and  D< 

o  •  -  p  -  D(1  +  cos  24)  -  C(24  ♦  sin  24) 

0#  -  -  p  -  D(1  -  cos  24)  -  C(24  -  sin  24)  (A-7) 

T  •  D  sin  24  +  c(l  -  CO*  24) 


It  is  now  convenient  to  express  the  «r  ,  o#  and  t  in  terms  of  principal  stresses 
,  e2  and  of  the  angle  7  -  7(4 )  which  the  direction  of  makes  with  the  radial, 
(Pig.  A-l).  Further,  the  stresses  and  «g  may  be  expressed  by  two  other  functions 
of  4,  namely  a  s  «(♦)  and  £  ■  £(4), 


-  -  («  +  •) 


*2  * 


-  £o 


(A-8) 


Using  the  conventional  expression  for  the  stresses  in  terms  of  the  principal 
stresses,  one  finds 

or  ■  cos^y  +  Og  sin27  »  -  [(l-£)e  ♦  s]  cos2y  -  £* 

0^  -  s^  sln2y  +  «gco827  -  -  ((l-£)o  ♦  s)  sin27  -  £s  (A-9) 

T  ■  («^  -  ag)sin  7  cos7  •  -  t(l-£)e  ♦  *1  sin  7  cos  7 
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Elimination  of  9  from  these  equations  gives  two  relations  between  ,  0^  , 
t,  k  and  7 


®r  -  «♦ 

t  -  ---g—  tan  27  -  0 
t  -  ■— -g~*  (—)  sin  27  +  sin  27  -  0  . 


(A-10) 


a.  SOLUTIONS  FOR  WEDGE  SHAPED  REGIONS WHEN  /(A)  -  +  £  . 


Applying  these  relations  at  the  boundary  ♦  •  0,  using  the  values 

£(0)  -  k  ,  7(A)  -  +  J  (A-ll) 

one  substitutes  the  values  of  or  «  0^  and  t  for  ♦  ■  A  from  Eqs.  (7 ),  and 
obtains  two  simultaneous  equations  for  the  constants  C  and  D: 


C  sin  20  ♦  D  cos  20  ■  0 
C  [(1  -  cos  20)  +  20  ♦  D  [sin  20  ± 


*[»&♦&] 


(A-12) 


The  upper  and  lower  sign  is  to  be  used  in  this  and  subsequent  equations  as  in 
Eq.  (11).  The  values  of  the  constant  C  and  D  are 


where 


C»  N  cos  20 
D  -  -  N  sin  20 


N  -  + 


„  1*k  *  k» 
P  IS  + 


1+k 


1-cos  20  ♦  (sin  20  -  20  cos  20) 


(A-13) 


(A-1A) 


It  is  noted  that  the  denominator  of  this  expression  may  vanish  for  certain 
combinations  of  0  and  k,  if  the  lover  alga  applies.  In  such  cases,  no 
solutions  exist. 
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It  is  nov  necessary  to  investigate  if  any  of  the  solutions  obtained  satisfy 

the  Coulomb  rule  Eq.  (3).  It  is  easy  to  see  that  for  the  special  values 

3  ■  £  and  3  -  ^  x,  Eq.  (3)  is  satisfied  everywhere.  In  these  cases  C  •  0  , 

and  the  state  of  stress  defined  by  Eqs.  (7)  is  uniform,  o  and  a  being 

x  y 

constant,  while  t  vanishes.  Having  prescribed  that  the  slip  condition  (2) 
is  satisfied  at  ♦  -  3,  it  is  obvious  that  Eq.  (3)  is  satisfied  everywhere. 
This  gives  the  following  simple  solutions 


3  ■  £  *  7(3)  m  +  % 
3-$,  7(3) 


x 

H 


-  p  ,  «x  - 


k(p-0*  t  -  0 


3  •  *  7(3)  •  ♦  £ 
3  -  £  ,  7(3)  m  ~  \ 


(A-15) 


-  p  > 


-  (|  ♦  •)#  t  -  0 


The  discussion  of  the  solutions  for  other  angles  3  becomes  much  more  Involved. 

As  a  first  step,  the  angle  7  can  be  determined  from  Eq.  (10)  by  substitution 
of  Eqs.  (7)  and  (13 )s 

cot  {&7)  •  •  <***> 
The  fact  that  the  denominator  vanishes  for  ♦  •  0  and  ♦  ■  23  -  a  (if  3  >  g) 
presents  no  difficulty.  One  can  now  determine  the  derivative  7'  "  ^  .  As 
will  be  seen  shortly,  the  behavior  of  the  quantity  1+7'  is  of  Interest,  and 
one  finds 


1+7' 


cos  23  [cos  23  -  cos  2(3-01 


1  +  cos  (23)  -  2  cos  (23)  cos  2(3-#) 
Applying  this  equation  for  the  boundary  ♦  ■  3  one  finds 


1  -  cos  23 


(A-17) 


(A-18) 
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In  the  range  0  <  0  <  a  the  right  hand  aide  haa  the  algn  of  (-  coa  20)  and  one 
find8  the  inequalitlea 


1+7' (0)  <  0 


(A-19) 


J<0<^  ,  1 +r(0)  >0  .  (A-20) 

d£ 

Aa  a  next  step,  an  expression  for  the  derivative  £'  c  ^  will  be  obtained. 

For  thiB  purpose,  Eqs.  (9)  are  substituted  into  Eqs.  (4)  giving  after  manipulations 

(1+y*  )(e(l-£)  ♦  a]  ain  y  +  £•'  coa  y  ♦  *£'  coa  y  «  0  (A-21) 

(1+7' )(«(!-£)  +  a]  cos  7  ain  y  -  0  .  (A-22) 


Excluding  locations  vhere  sin  2y  -  0,  one  can  eliminate  •'  and  obtain  a  relation 

(A-23) 


£•  .  (u£)  [;lV>i  [co.  2r  -  ii] 


Applying  this  result  to  the  determination  of  the  algn  of  £'  at  ♦  •  0,  one  has, 

In  addition  to  Eqa.  (19,20),  the  following  Information 

£(P)  •  *  t  7(0)  »  ♦  £ 

l+£  >0  ,  s(l-E)  ♦  a  >  0  ,  a  >  0  (A-24) 


coa  2y  - 


k! 

i+£ 


l-k 

1+k 


<  0 


Using  Eqs.  (19),  (20),  (23)  and  (24)  one  finds  that  £'  <  0  provided i 
7(P)  •  £  and  {•  <  0  < 


or 


(A— 2J5) 


7(0)  -  -  J  *nd  0  <  0  <  j|  or  <  0  <  «  . 
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For  the  complementary  range*  of  0  one  finds  £'  >  0.  The  Important  point  is 
that  the  material  at  the  boundary  ♦  •  0  is  Just  at  the  verge  of  the  slip, 
k  ■  k,  such  that  the  condition  £'  <  0  must  be  satisfied,  otherwise  the  Coulomb 
rule  is  violated  at  small  distances  from  ♦  •  0.  This  restricts  the  ranges  to 
be  studied  further  to  those  defined  in  Eq.  (25 )• 

In  the  geometries  defined  by  Eq.  (25),  the  Coulomb  rule  is  satisfied  near 
♦  •  0j  to  investigate  other  locations,  Eq.  (16)  is  used  to  obtain  the  range  of 
7  as  ♦  varies  from  0  to  0.  Equation  (A-l6)  is  simple  enough  to  drav  general, 
qualitative  conclusions.  The  three  cases  in  Eq.  (25)  separate  into  a  total  of 
four  situations,  listed  in  the  following  table  depending  on  the  prescribed  value 
of  7(0)  ■  +  jj  and  on  the  wedge  angle  0.  The  function  >(♦)  decreases  In  all 
cases  monotonlcaliy  from  7(0)  to  7(0): 


(A-26) 


To  proceed,  a  general  expression  for  £(♦)  is  required,  and  may  be  obtained  from 
Eqs.  (10).  On*  could  solve  the  second  of  these  equations  for  £,  but  the  result 
would  break  down  whenever  sin  27  -  0,  l.e.,  according  to  Iq.  (26)  on  the  surface 
♦  m  0.  To  circumvent  this  difficulty,  one  forms  the  difference  of  the  two  equa¬ 
tions  (10)  and  finds  the  alternative  relation 


l-£  *r  *  £s 

±-£  -E - L.  cos  2r  ♦  ==-  cos  27-O 

l+£  2  l+£ 
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(A-27) 


Using  Eqs.  (7)  and  (13)*  one  obtains 


P+N  ’  *ln  *  2*  008  2p] 
p—N  ["--co^r1  +  Bin  2*  - 2* COB  2(J] 


(A-28) 


This  result  can  be  used  to  check  on  the  values  of  £(0).  Taking  tha  value  of  7(0) 
from  Eq.  (26),  and  considering  values  0  permitted  by  Eq.  (2$)  one  finds: 

If  0  <  0  <  ^  *orif^<0<^p,'  the  value  of  the  cosine  is 

cos  [2y(0)]  »  +  1,  and 


£<°>  -  p— -SH^.lo  <*-*»> 

If*  however,  £  <  0  <  |  ,  or  if  ^  <  0  <  *  *  then  cos  (27(0)]  ■  -  1,  end 

£(0)  -  (A- 30) 


Before  evaluating  Eqs.  (29*30)  knowledge  of  the  permissible  range  of  £  is  required. 
From  the  definition,  Eq.  (8),  for  the  principal  stresses  it  is.  clear  that  £  >  k 
is  a  sufficient  condition  for  the  satisfaction  of  the  Coulomb  rule  if  0^  is  the 
major  stress.  But  there  will  be  locations  where  e^  is  the  major  stress*  and  in 
such  locations  £  must  not  exceed  an  upper  limit*  which  can  be  easily  computed. 

If  s  •  0,  the  upper  limit  is  simply  i  *  but  in  general 

*  s  £  <  s  ♦  ^  <*-»> 

figure  (A-2)  gives  the  typical  appearance  of  the  plot  of  the  four  branches  of  £(0). 
Regardless  of  the  values  of  s*  the  value  of  £(0)  is  always  at  the  limit  of  slip 
if  0  ■  ^  or  *  because  the  solution  Eq.  (1$),  (a  constant  state  of  stress)  applies 
there.  For  0  <  jj  or  0  >  (Fig.  A-2)  shows  that  the  Coulomb  rule  is  always  vio¬ 


lated*  and  no  solutions  exist  for  such  angles. 


In  the  range  ^  <  0  <  ^  the  value  £(0)  has  a  maximum  for  0  ■  |  .  To 
verify  that  Eq.  (31)  la  satisfied  at  this  point,  Eqs.  (13)  are  substituted  in 
Eqs.  (7)  and  one  finds  for  ♦  ■  0,  of  ■  •  p  ,  t  •  0,  *  state  of  stress  vhlch 

certainly  satisfies  the  Coulomb  rule  and  therefore  Bq.  (31).  If  s  ■  0,  the 
peak  value  of  £  at  0  »  ^  is  unity. 


Ve  have  nov  arrived  at  the  conclusion  that  In  addition  to  Eqs.  (15) 
suitable  solutions  may  exist  only  If  y(0)a  +  ^  and  ^  <  0  <  .  It  Is  still  not 

known  If  the  solutions  in  this  range  satisfy  the  Coulomb  rule  expressed  by  Bq.  (31) 
everywhere  or  not.  To  confirm  that  (31)  Is  satisfied  it  Is  necessary  to  consider 
the  situations  for  0  ^  |  separately.  Tor  simplicity,  the  argument  will  be  demon¬ 
strated  for  s  ■  0  only. 

For  the  case  ^  <  0  <  g  the  value  of  y  ranges  from  ^  to  j  (for  ♦  •  0),  such  that, 
with  the  exclusion  of  ♦  -  0, 


sin  2y  >  0  and  cos  2y  <  0  . 


In  the  above  range  for  0  the  sign  of  1  +  7'  given  by  Eq.  (17)  does  not  change, 
and  for  ♦  jfc  0 

1  +  7'  >  0  . 


The  ratio 


1  +  7* 
sin 


fr 


>  0 


is  therefore  positive;  this  statement  holds  even  In  the  limit  ♦  •  0. 
Equation  (23)  gives  the  sign  of  £'  for  s  ■  Ot 

sign  (£• )  -  sign  |^l-£^  ^cos  2y  -  • 


(A- 38) 


It  is  known  already  from  Fig.  (A-2)  that  1  >  £(0)  >  £(0)  ?  k,  and  the  above 
equation  indicates  that,  beginning  at  ♦  ■  0,  the  value  of  £'  Is  negatlvs. 
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As  ♦  increase*  this  could  only  change  if  k  became  equal  to  unity,  which  is 
impossible  as  k(0)  <  1.  The  value  of  £  will  therefore,  decrease  smoothly  frost 
♦  ■  0  to  ♦  «  fl,  such  that  the  Coulomb  rule  for  ^  <  0  <  |  is  satisfied  everywhere. 

For  values  |  <  0  <  the  proof  is  a  little  more  complicated,  because 
£(♦)  increases  from  ♦  ■  0  to  a  maximum  value,  and  subsequently  decreases  to 
i(P)  •  k. 


According  to  Eq.  (26),  7  ranges  from  y(0)  -  *  to  7(0)  -  ^  .  Equation  (16) 
yields  easily  the  locations  where  7  ■  ^  or  ^  ,  respectively: 

7(P  *  f)  -  and  7 ^  -  *)  “  \  • 

Using  this  knowledge  one  can  show  again  that 

— -Zl  >  o 
sin  27 

for  the  entire  range  0  <  ♦  <  0  ,  and  the  sign  of  £'  is  therefore  again  given  by 
Eq,  (32).  Substituting  the  appropriate  values  of  7  at  ♦  •  0  ,  0  and  noting  that 
in  these  locations  k  <  1  ,  one  finds 

£' (0)  >  0  ,  £'(0)  <  0 

Indicating  a  maximum  of  £  in  the  interior.  To  locate  the  maximum,  examine  the 
value  of£at  +  -0>|,  where  7  -  and  cos  27  •  0  ,  sin  27  »  »  1. 

In  general,  the  value  of  £(♦ )  is  given  by  Eq.  (28),  but  this  equation 
breaks  down  if  cos  27  -  0.  However,  using  Eq.  (l6)  one  finds 


sin  2(p-»)  ^  cos  20  -  cos  2(0-») 
cos  27  "  sin  £7 


(A-33) 


Substituting  this  expression  in  Eq.  (28),  one  obtains  a  valid  expression  for 
%  i.C 

♦  ■  0  -  5  .  Forming  the  combination  — —  ,  one  obtains  (for  s  ■  0) 

*  i+£ 
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l+k 


>  0  . 


(A-3M 


(1  ♦  eo*  20) 

m 

x  1  -  cos  20  -  x  T-r  cox  20 
2 

The  Inequality  lx  easily  proved,  becauxe  in  the  range  conxldered  ve  have 
0  >  cox  20  >  -  1. 

The  above  Inequality  igqpliex  £  <  1,  and  Eq.  (32),  with  cox  2 y  ■  0,  givea 

k'l  «  <0 

!♦  -  0  -  5 

indicating  that  the  maximum  value  of  k  occura  between  ♦  ■  0  and  ♦  ■  0  -  j  ,  In 
this  range  the  angle  y  ia  Halted  by^p<y<xor0<  cox  2  y  <  1. 

For  the  maximum  value  of  £  the  derivative  (32)  must  vanish.  Approaching  from 
♦  -  0,  having  started  with  £  <  1,  and  noting  the  above  Inequality  on  cox  2y,  it 
lx  clear  that  the  xecond  factor  in  Eq.  (32)  will  be  the  one  to  vanish  when 
£  -  max.  £  ;  or 

_  1  -  aax  £ 

cox  zy  *  . . . .  •  0  . 

1  ♦  max  £ 

In  this  equation  cos  2y  >  0,  such  that  aax  £  <  1  completing  the  proof. 

In  the  final  stage,  only  the  case  a  ■  0  has  bean  considered.  However,  the 
values  of  £  being  a  continuous  function  of  s,  one  reasons  that  the  conclusions 
must  also  hold  for  small  values  of  the  ratio  j  (small  vs.  unity)/  This  region 

suffices  for  the  present  purpose,  as  values  of  p  comparable  to  s  are  not  of 

•  ) 

Interest  anyway  . 

*)  It  is  demonstrated  in  subsection  b  that  elastic  solutions  exist  up  to  some 
value  £  >  1,  the  critical  ratio  being  a  function  of  k.  Xf  solutions  for 
wedges  without  slip,  satisfying  the  Coulomb  rule  everywhere  did  exist  for  all 
values  of  ^  <  0  <  ^  ,  even  for  values  where  j  is  not  small  vs.  unity,  they 
would  ultimately  lead  to  solutions  (for  the  half-space)  of  the  steady-state 
problem  with  slip,  for  values  p  where  elastic  solutions  also  sxlst.  Xn  such 
a  case  solutions  with  slip  would  not  be  expected  to  have  physical  significance, 
even  if  they  do  exist. 
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Suanarltlng,  it  bu  been  found  that  for  y(0)  -  £  solution*  exist  only 
for  wedges  J  <  0  <  %  ,  vhlls  for  y(0)  -  -  {  solutions  (Iq.  15)  sxlst  only 
for  0  •  or  p  -  « 

For  the  ease  y(0)  ■  ^  ,  the  following  expreeelon  for  the  value  of  the 
stress  pereaeter  e  on  the  Inclined  face  ♦  -  0  can  be  obtained! 

2p(l  -  cos  20)  -  s(l  -  cos  20  •  20  cos  20  ♦  sin  20) 

•(l)  ’  cuk)(i  ^csrsgTTirEKHrSg :  cos  20 )  ^  (A-35) 
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b.  80LUTI0N8  WITHOUT  SLIP  TOR  THE  HALT-SPACE. 

In  order  to  determine  up  to  what  level  of  tbe  load  p,  the  elastic  stresses 
satisfy  the  Coulomb  rule  for  a  half-space,  the  left  side  of  vhlch  Is  loadsd, 

(Fig,  A-3),  the  solution  Is  obtained  from  Bq.  (6).  Surprisingly,  the  four 
boundary  conditions,  e^(0)  ■  0,  s^(s)  -  -  p,t(0)  •  t(k)  ■  0  permit  the  deter¬ 
mination  of  only  three  of  the  four  constants.  The  last  one,  D, remains  open, 

such  that  the  stresses  become,  with  6  -  -  -  D 

P 

—  -  -  ♦  -  ^  sin  2#  ♦  B(1  ♦  cos  2*) 

-  -  ♦  ♦  i  sin  2#  ♦  D(1  -  cos  24)  (A-36) 

P  • 

-  ^(1  -  cos  2*)  *  C  sin  2*  . 

A  restriction  on  the  constant  9  follows  Immediately  by  considering  the  stresses 
at  the  free  surface  ♦  -  Oj  e^  •  x  •  0  and  ey  ■  .  Requiring  one  principal 

etreas  to  be  compressive,  the  Coulomb  rule  gives 

0  >  *  *  -  g  •  (A- 37) 

The  free  constant  is  physically  equivalent  to  an  addltlonsd  uniform  horizontal 
state  of  stress  0  -  ,  s  ■  t  •  0.  The  lack  of  uniqueness  of  the  solution 

Is  a  consequence  of  the  lack  of  a  prescribed  state  of  stress  at  Infinity. 


To  determine  if  slip  occurs  at  any  point,  tbe  value  of 


£  . 


•l  +  * 


is  required.  Using  the  usual  expressions  for  the  principal  stresses  one  obtains 


£  . 


D 


•  -s  +  TfS'-f 


(A-36) 
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where 


A  -  Bit?  ♦  ♦  f?  -  2S  ala  ♦  ooa  ♦  (A-39) 

If  the  value  of  IT  la  given,  one  could  find  the  minimum  value  of  E  for  any 
ratio  £  ,  and  the  aolutlon  would  aatlafy  the  Coulomb  rule  when  k  <  aln  E  . 

A  .change  of  5  will,  of  course,  effect  ein  E  ;  it  ia  of  intereet  to  know  the 
•oat  beneficial  value  of  IT  which  raleee  aln  E  to  ita  highest  value,  because  It 
defines  values  £  above  which  elastic  solutions  without  slip  cannot  exist* 


To  find  extremal  values  of  £  with  respect  to  ♦  and  5,  the  somewhat  simpler 
expression  1  t  E  is  formed, 


„  2(B-#)  ♦  H 

1+E  - - 2-- 

and  differentiation  with  respect  to  ♦  and  B  gives  two  equations 

[«»*>  ♦  *]  ['  1  -  ft]  •  *  «»  *  *  -  1?*  *> 


(A-ko) 


(A-kl) 


Addition  of  these  equations  indicates  that  either 


2(B-*)  +  SL  .  o 

or 


(A-ka) 


dA  .  dA 
5?  +  3B’  0  * 


(A-kj) 


If  lq.  (k2)  applies,  the  right  hand  side  of  (kl)  must  also  vanish,  giving 


♦  -  B  -  /a1,  and  according  to  lq.  (3d),  E  »  0.  This  situation  is  trivial. 
The  alternative,  lq.  (k3),  becomes  after  substitution 


k  B  sin2  ♦  m  0 

yielding  the  conditions 

D-0,  or  ♦  ■  0,  or  ♦  ■  *  . 


(A-kk) 

(A*kJ) 
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The  last  two  possibilities,  ♦  ■  0  or  *,  lead  again  to  trivial  cases. 


This  leaves  the  root  5  ■  0.  Substitution  of  this  value  in  either  of 
Eqs.  (Ul)  gives  a  relation  between  j  and  ♦  : 


£ 


*(1  -  cos  ♦) 
2(sln  4  -  ♦  cos  4) 


(A  -*6) 


It  can  be  shown  that  the  solution  mini  mites  £  with  respect  to  ♦  ,  and  that 
fi  •  0  gives  the  highest  values  of  the  minimum  as  desired.  Computing  min  £ 
and  equating  it  to  the  value  k  of  the  material  to  be  considered,  one  finds 


k  -  sin  £ 


1  -  cos  ♦ 

1  +  cos  ♦  * 


(A-Vf  ) 


This  is  a  parametric  representation  of  £  as  function  of  k.  It  defines  the 
values  of  p  shove  which  no  elastic  solution  without  slip  exist,  figure  (A-U) 
shows  a  plot  of  £  versus  k.  It  is  noted  that  the  location  ♦  at  which  the 
minimum  occurs  varies  with  k;  for  small  values  of  k  the  angle  ♦  is  also  quite 
small;  ♦  increases  with  k,  and  reaches  the  value  ♦  •  |  for  k  ■  1. 
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APPENDIX  B 


ON  THE  DISINTEGRATION  OF  BODIES  OP  MATERIALS  GOVERNED  BY  THE  COULOMB  RULE. 

In  the  materials  considered  the  premissible  states  of  stress  are  limited 
by  the  inequality  (2-3) 

-  *2  >  k«. 

If  one  subjects  a  body  to  loads  where  this  condition  cannot  be  satisfied, 
the  body  being  unable  to  support  the  loads  must  disintegrate  entirely  or  at 
least  in  part.  To  establish  the  manner  in  which  the  disintegration  occurs 
requires  a  physical  description  of  the  behavior  of  the  elements  of  the  body 
beyond  the  mathematical  statements  in  Section  I.  In  the  following  a  simple 
situation  will  be  considered  in  which  the  material  can  be  represented  by  a 
mechanical  model  which  has  properties  in  agreement  with  relations  Eq.  (2-2  to 
2-8)  when  the  cohesion  s  vanishes.  (An  alternative  model  with  cohesion  can  also 
be  formed). 

Consider  the  plane  structure  covering  a  rectangular  area,  (Pig.  B1 )j  the 
structure  is  built  of  Identical  square  elements  and  carries  loses  Pj  ,  in  one 
dlrectlon,and  Pg  in  the  other.  Each  element  consists  of  four  masses  connected 
by  lnextenslble  linkage  bars  under  45  degrees  to  the  vertical,  (Fig.  B2).  Let 
the  elements  be  small  and  of  width  2a  and  height  2a.  The  hinges  are  assumed  to 
have  friction,  the  friction  moments  being  proportional  to  the  load  transferred 
to  the  hinge  by  the  bar.  If  the  unit  is  to  be  in  equilibrium  under  forces 

>  Pg  acting  on  the  element,  (Pig.  B2),  the  forces  must  satisfy  the  Inequality 

Px  >  P2  >  kPx  (B-l) 

If  however 

P2  -  kPx  (B-2) 

the  unit  can  be  in  steady  motion  or  at  rest.  The  coefficient  k  can  be  obtained 
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by  equilibrium  consideration*  from  the  friction  in  the  hinges,  where  by  reason 
of  synmetry  the  forces  in  one  link  only  need  be  considered. 

The  friction  in  a  hinge  may  be  defined  by  the  maximum  normal  distance  f 
at  which  the  resultant  force  8  in  the  linkage  bar  will  bypass  the  center  of  the 
hinge.  This  applies  to  the  hinges  at  both  ends  of  the  link  and  due  to  the  type 
of  linkage  motion,  the  resultant  must  pass  through  the  mid-point  of  the  linkage, 
(fig.  B-3).  The  resultant  will  therefor*  ask*  an  angle  a  with  the  link,  such  that 

*  -  <B-3> 

Later,  it  will  be  convenient  to  know  the  connection  between  a  and  k.  Equilibrium 
at  the  four  hinges  requires  the  following  two  relations  between  the  compressive 
resultant  6  and  the  forces  and  Pg  -  kP^  ,  acting  on  the  eleswnt  in  the  limiting 

state, 

P^  ■  28  cos  (J  -  a) 

P2  -  kPx  -  28  sin  ({  -  a) 
or 

k  •  tan(^  -  a)  .  (B-4) 

The  changes  in  vertical  and  horizontal  "strain"  ~  ,  (Pig.  B-2),  will  be  of 
equal  magnitude  in  similarity  with  Bq.  (2-8).  It  is  further  aaeuaed  that  the 
angular  motion  of  the  links  is  restricted  to  a  small  angle,  such  that  the 
motion  a  cannot  exceed  a  set  value  a  .  When  this  situation  is  reached,  the 
unit  will  be  deemed  to  become  rigid. 

The  element  described  behaves  like  the  material  considered  under  A  in 
Section  2,  provided  the  vertical  axis  of  the  element,  (Pig.  B-2),  is  in  the 
direction  of  the  major  principal  stress.  The  only  difference  lies  la  the 
limitation  on  the  strain  resulting  from  the  condition  a  <  a  . 
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A  number  (2N)  of  such  element*  arranged  in  a  horizontal  line,  (Fig.  B-h), 
will  be  in  incipient  motion  if  each  element  la  vertically  loaded  by  the  force 
P,  while  horizontal  force*  kP  are  applied  at  the  enda  and  tranamltted  from 
element  to  element.  The  arrangement  almulatea  a  rectangular  region  of  a  Coulomb 
material  at  the  verge  of  alip.  Let  ua  determine  the  reaponae  of  the  model  if 
the  force*  kP  on  the  aide*  are  auddenly  removed.  The  behavior  of  the  modal  ought 
to  give  an  inaight  on  what  will  happen  in  caae  of  a  granular  material. 

It  baa  already  been  atated  that  the  linkage  bar*  are  lnextenslble;  for 
aimpliclty  it  ia  further  aaaumed  that  they  are  massless,  while  each  of  the  four 
maeeea  at  the  hinge a  baa  one  quarter  of  the  total  aaaa  M  of  each  unit.  Also, 
for  alnqpllcity,  it  i*  aaaumed  that  rotation  of  the  maaaea  la  prevented  (e.g.  on 
top  and  bottom  by  the  atlffneaa  of  platea  in  a  teating  machine  by  me ana  of  which 
the  loada  P  are  applied,  while  the  other  maaaea  will  not  rotate  by  reaaon  of 
eymmetry). 

Conaider  the  motion  of  the  1-th  element,  (Fig.  B-5 ).  It  la  acted  on  by 
the  external  force*  P  on  top  and  bottom;  by  aa  yet  unknown  horizontal  fcrcaa 
and  P^  from  the  two  adjoining  element*;  in  addition,  at  each  end  of  each 
linkage  bar  there  will  be  friction  momenta  acting  on  theme  bare,  f8^  and  f8g  , 
where  8^  and  8^  are  the  reaultanta  transferred  by  the  left  and  right  linkage, 
respectively.  At  the  top  and  bottom  binge*,  the  friction  moments  do  not  balance, 
the  difference  being  supplied  by  the  external  support  preventing  rotation. 

The  notion  of  the  element  is  fully  described  by  two  generalized  coordinates! 
the  horizontal  displacement  x^  of  the  center  0^  of  the  unit,  and  the  relative 
motions  ,  horizontal  and  vertical,  respectively,  of  the  mass  points  with 

respect  to  the  center  0^  .  For  reasons  of  eymmetry  in  the  arrangement,  (Fig.  B-b), 
no  vertical  notion  of  the  center  0^  need  be  considered.  For  small  motions,  the 
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rotation  of  the  linkages  is  |  .  To  write  Lagrange's  equations;  one  obtains  tbs 
kinetic  energy 


2T  .  N  (xj  +  «J) 


(B-5) 


while  the  potential  energy  vanishes.  The  generalised  forces  Q  and  Q  due  to  the 

xi  fi 


physical  forces  shown  in  (Fig.  B-5)  ore 


V  -  Pi-1  -  ?i 


2P 


hr 


(B-6) 


Pi-1  -  Pi  •  T  <*L  ♦  V 


Noting  Kq.  (3),  the  equations  of  notion  become 


**!  "  \  *  Pi-1  ’  Pi 


(B-7) 


¥k\  -  -  2P  -  P1.1  -  Pt  -  2  (F sin  a  (SL  ♦  Bg).  (B-8) 


The  second  equation  still  contains  the  resultants  8^  and  8^  in  the  linkage 
bars,  which  can  be  expressed  in  terns  of  3^  and  .  The  forces  acting  on 
the  two  lateral  aasses  are  shown  in  (Pig.  B-6);  their  accelerations  are 
x^  +  ,  respectively;  such  that 

5  (*1  ’  «i)  "  pi.i  "  2  Bln  ({  *  *) 


*  ♦  i’j)  ■  2  sin  "  *)  Bg  “  P^ 


or 

2  sin  -  aHBg  ♦  8^)  -  |  «t  ♦  (P4  ♦  P^).  (B-9) 

Substitution  into  Iq.  (8)  gives 


»  -  t  (»1  ♦  ?l-l> 


(B-10) 
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where  the  relation 


aln  a _  1-k 

$  ■ln(*  -  a)  *  2k  (B-ll) 

TT 

was  obtained  from  Eq.  (U). 

Eq.  (10)  vaa  derived  on  tha  i»cit  understanding  that  g^  >  0,  indicating 
that  the  element  deforma  in  the  manner  shown  in  (Fig.  B-5).  An  alternative 
relation  valid  if  <0  could  easily  be  derived,  but  is  not  required  for  the 
present  purpose.  Further,  it  is  possible  that  the  external  and  inertial  forces 
are  such  that  tha  friction  is  sufficient  to  prevent  any  deformation  of  the  unit, 
such  that  ■  0.  If,  in  this  case  the  unit  is  not  subject  to  acceleration, 

■  0,  the  forces  and  vill  be  equal, 

P  ■  P 

1  1-1  (B-12) 

and  must  satisfy  the  Coulomb-inequality 

i  P  >  P.  >  kP  (B-13) 

k  1 

Equivalent  relations  for  x^  i  0  could  be  derived,  but  vill  not  be  required  for 
the  present  purpose. 

The  differential  equations  (7)  and  (10)  and  the  relations  (12,  13)  permit 
determination  of  the  response  of  the  group  of  elements  shown  in  (Fig.  B-4)  when 
the  lateral  loads  kP  are  suddenly  removed  at  a  time  t  ■  0.  Considering  x^,  Sj 
and  P^  as  functions  of  t,  the  forces  Pg  and  P_g  must  vanish  for  t  >  Oi 

P  l  P  -  0  (B-ll*) 

H  -N 


Further,  when  adjoining  elements  are  in  contact,  the  geometric  relation 


!♦! 


-  x. 


«!♦!  *  ‘i 


(B-15) 
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must  b«  satisfied.  At  t  ■  0,  Eq.  (15)  will  apply  for  all  values  X-l  >  1  »  -I, 
but  thereafter  separation  night  occur  If  *  0,  or  due  to  the  limitation  on 
the  notion  of  the  linkages,  <  f.  These  possibilities,  and  the  question 
whether  the  strain  rate  is  positive,  negative  or  tero,  nust  be  carefully 
considered. 


To  start,  the  sinple  ease  of  two  elements,  B  ■  1,  will  be  solved,  (Fig.  B»T)» 
where  for  reasons  of  symmetry  only  unit  (l)  need  be  studied.  Subject  to 
later  check,  it  is  assumed  that  >  0,  such  that  Eqs.  (7)  and  (10)  apply. 


•  P- 


1+k 

2k 


•  2  P 


(B-16) 


Eq.  (15)  becomes 


and  elimination  of 


and 


gives  the  simple  relation 


(B-17) 


(B-l8) 


P  is  positive  as  required.  Further,  using  the  initial  conditions 
o 

*  •  ■  a  ■  *  •  0,  one  finds 

*1  >  °*  *i  >  °*  *1  *  °*  *i  *  ° 

For  sufficiently  small  values  of  t  all  conditions  for  the  validity  of 
Eqs.  (16)  are  therefore  satisfied.  The  centroids  0^  of  the  elements  move  away 
from  each  other,  while  the  "strain  rate"  is  positive.  However,  at  an 
instant  t  ■  ?,  the  strain  reaches  the  limit  s ,  and  the  elements  freese. 


6l 


At  this  instant,  the  centroids  of  both  elements  still  have  outward  velocities, 


x  > 
1 


such  that  the  elements  will  separate  and  fly  apart. 


The  value  of  the  pressure  between  the  two  elements  prior  to  separation 
is  given  by  Eq.  (l8).  It  is  noted  that  the  value  PQ,  regardless  of  the  exact 
value  of  k,  is  always  slightly  larger  than  the  value  kP  required  to  prevent 
<lip  under  static  conditions.  This  means  that  the  solution  for  a  larger  number 
of  elements,  N  >  1,  can  be  stated  without  further  analysis.  For  t  <  t,  the 
outside  elements  N,  -N,  in  Fig.  (B— U )  behave  exactly  as  the  element  Just 
analyzed  in  the  case  N  •  1.  The  interior  elements,  N-l  to  -(Hoi)  do  not  move 
or  change  shape,  because  the  force,  Eq.  (18),  exerted  by  the  outside  elements 
does  not  produce  slip.  At  t  ■  F,  the  outside  elements  separate,  and  a  new 
situation  arises.  The  elements  (N-l)  and  -(N-l)  begin  to  move,  while  no  motion 
occurs  in  the  other  elements,  until  separation  occurs  again,  etc. 


We  find,  therefore,  that  in  case  of  a  large  number  of  elements,  the  dis¬ 
integration  of  the  model  structure  occurs  in  a  layer  of  elements  at  the  boundary. 
The  accelerations  of  the  elements  in  the  boundary  layer  produce  a  reactive  pressure 
which  produces  in  the  Interior  a  state  of  stress  satisfying  the  Coulomb-rule? 


In  the  above  derivations,  linearity  was  assumed,  and  enforced  by  permitting 
only  a  finite  small  deformation  7  of  the  elements.  If  the  restriction  had  been 
dropped,  the  element,  after  experiencing  large  strains,  could  still  not  exceed  a 
certain  finite  strain  vhere  the  linkage  is  fully  extended|  separation  would 
therefore  occur,  and  is  not  Just  due  to  the  assumption  <  c. 

Visualizing  the  mechanism  of  acceleration  and  subsequent  separation  in  a 
granular  medium,  one  can  see, (Fig.  B-8),  how  a  grain  3  might  be  expelled  by 
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Appropriate  wedge-like  surfaces  of  two  other  grains,  1  and  2,  if  the  latter 
are  pressed  together. 

In  analogy  to  the  nodal,  ve  expect  that  on  a  free  surface  in  a  Coulomb- 
type  granular  material,  disintegration  will  occur  if  the  problem  has  no 
solution  for  which  the  pressure  acting  on  the  surface  vanishes.  The  proceee 
of  disintegration  and  acceleration  will  occur  in  a  thin  boundary  layer,  and 
will  restore  a  pressure  sufficient  to  satisfy  the  Coulomb  inequality  below 
the  surface. 

One  might  be  tempted  to  critise  the  use  of  the  mechanical  modal  because 
it  cannot  represent  a  Coulomb  material  if  the  principal  stresses  should  change 
direction  during  loading.  However,  this  criticism  Is  not  pertinent  because 
the  process  of  disintegration  occurs  in  the  immediate  vicinity  of  a  surface 
where  such  directions  necessarily  eannot  change. 

The  model  could  be  adapted  to  permit  representation  of  a  material 
with  dllatancy,  or  with  cohesion.  In  the  former  ease,  linkages  making  angles 
other  than  90°  would  be  used.  A  model  for  the  seoond  oase  is  obtained  by 
specifying  the  friction  in  the  hinges  to  be  a  constant  plus  f8.  Such  a  model 
could  also  be  given  tensile  strength  by  providing  attraction  between  the 
elements.  The  basic  conclusions  are  not  af footed  by  the  use  of  the 
above  modifl cations. 
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In  Appendix  A:  arbitrary  constants  of  differential 
equation. 

Acceleration  in  radial,  x-  and  y-  direction,  respectively. 

In  Appendix  B:  one  half  of  height  of  elements. 

Velocities  of  compressive  and  shear  waves,  respectively, 
in  an  elastic  material. 

In  Section  1:  a  constant  in  Eq.  (1-1). 

Basis  of  natural  logarithms. 

Positive  quantity,  defining  the  principal  strain  rates. 

Young's  modulus. 

In  Section  1:  yield  function,  plastic  potential,  respectively. 
In  Section  3:  any  function. 

In  Appendix  B:  distance  shown  in  (Fig.  B-3). 

In  Section  1;  an  arbitrary  function. 

Modulus  of  rigidity. 

Subscripts,  i,J  «  1,2,  ... 

Invariants  of  strain  deviator. 

Constant  defining  the  stress  at  which  slip  occurs  in 
Coulomb  material,  Eq.  (2-6). 

In  Appendix  A;  a  parameter  expressing  the  principal 
stresses  in  a  non-slip  elastic  region,  Eq.  (A-8). 

Length  of  load  in  Fig.  3» 

Exponents  in  Eq.  (3~7)> 

In  Appendix  A:  expression  defined  by  Eq.  (A-14). 


*> 


Other  symbols  which  are  locally  used  are  defined  where  they  occur. 
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In  Appendix  B:  centroid*  of  elements. 

Zn  Appendix  B:  forces  acting  on  the  elements. 

Intensity  of  applied  step  pressure. 

Surface  pressure  at  ♦  ■  0  (ahead  of  applied  load). 

Limiting  value  of  p,  Eq.  (4*6). 

Maximum  static  pressure  according  to  Prandtl  [Ref.  2]. 
Limiting  value  of  p  below  which  solutions  without  slip  exist. 
Generalized  forces. 

Parameter  defining  cohesion. 

Zn  Appendix  B:  force  in  links. 

Time. 

Kinetic  energy. 

Components  of  velocity,  horizontal,  tangential,  radial 
and  vertical,  respectively. 

Velocity  of  applied  pressure  pulse. 

Cartesian  coordinates,  (Fig.  1). 

Zn  Appendix  B:  displacement  of  i-th  elements. 

Derivative  with  respect  to  t,  (,  respectively. 

Zn  Section  1:  Coefficient  in  Drucker  and  Prager's 
(Eq.  1-1),  also  Inclination  of  pressure  front,  (Fig.  4). 

Zn  8ection  2:  Angle  defining  position  of  slip  plane, 

(Fig.  5). 

Zn  Appendix  B:  angle  between  resultant  force  and  link, 

(Fig.  B-3). 

Opening  angles  of  wedge  shaped  regions  in  which  no  slip 
occurs. 

Angle,  cos  2y  -  • 
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Dirac's  Delta  Function. 

In  Appendix  A:  expression  defined  by  Eq.  (A-39)* 
Strain  and  strain  rate,  respectively,  with  appropriate 
subscripts, ■ j  ,  etc. 

In  Appendix  B:  Deformation  of  element. 

Principal  strain  rates  in  direction  of  principal 
stresses,  o^  ,  0g  ,  respectively. 

Angle  defining  direction  of  major  principal  stress, 
(Fig-  5). 

Location  dependent  quantity  vben  deriving  strain  rates 
from  potential  functions,  Eqs.  (1-2),  (1-3),  (2-9)* 
Poisson's  ratio. 

Non-dimensional  variable. 


p  Mass  density  of  material. 

e  Independent,  necessarily  positive  variable,  defining 

the  principal  stresses,  Eqs.  (3**2). 

a  ,  Og  Major  and  minor  principal  stresses,  respectively. 

°x  '  *y  '  °r  '  *♦  '  9iJ  stress  components. 
t  Shear  stress. 

♦  Position  angle  of  an  element,  (Fig.  6),  measured 


clockwise  from  horitontal. 

Note  re  Numbering  of  Equations. 

Throughout  this  report,  equations  are  numbered  by  hyphenated  numbers,  such  as 
Eq.  (2-7),  (3-18),  or  (B-lU).  The  first  number,  2,  3,  B,  respectively.  Indicates 
that  the  equation  occurs  in  Section  2,  3,  or  in  Appendix  B. 

When  an  equation  is  referred  to,  the  full  number,  say  (2-7),  ia  quoted  if  the 
reference  occurs  in  a  different  section.  However,  if  an  equation  is  mentioned 
in  the  section  in  which  it  originally  appears,  only  the  second  half  of  the  number 
is  quoted;  In  Section  2,  Eq.  (2-7)  would  therefore,  be  quoted  simply  as  Bq.  (7), 
etc. 
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